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CHAPTER I 
INTRODUCTION 


Among the various problems of electromagnetic scattering by arbi¬ 
trary bodies, only a few simple shapes have led to exact solutions. 

The diffraction of an electromagnetic plane wave by a conducting circu¬ 
lar disk long remained unsolved. The presence of an edge led to much 
greater complications than, for example, the sphere which was solved by 
Mie[l] very early in this century. 

First, some limiting cases were derived for the disk problem. 
Approximate solutions were obtained by Kirchhoff diffraction theory 
for disks of large electrical circumference kd where k is the wave 
number and d the diameter of the disk and by Lord Rayleigh[2] for 
disks of small electrical circumference. After the second world war, 
extensive work was directed at the disk problem, and two complete gen¬ 
eral solutions were developed. Meixner and Andrejewski[ 3 ] used Hertzian 
vector potentials to compute the scattered field. They derived the form 
of the potential from the appropriate boundary conditions, and the uni¬ 
queness of the solution was insured by the edge condition developed by 
Meixner[4]. He showed that the energy density of the total field in the 
vicinity of an edge must be integrable. This insists that the compon¬ 
ents of the fields vary at most as s"’'S where s is the distance to the 
observation point from the edge. Bouwkamp[5] showed that the tangential 
electric field is zero on the edge, that it approaches this value as 
s 1 ' 2 , and that the tangential magnetic field remains finite. Tangen¬ 
tial, in this context, means parallel to the edge. This solution has 
been checked by numerous computations made in recent years. Since the 
development by Hodge[ 6 ] of an efficient method of computing the spher¬ 
oidal eigenvalues, even more efficient computers have allowed extensive 
comparison between experimental data and computational results with very 
good agreement. 

Shortly after Meixner, Flammer[7] derived another solution using 
oblate spheroidal vector wave functions to expand the fields. As in the 
Mie solution for the sphere, the fields of the incoming plane wave were 
expanded in terms of the vector wave functions. The boundary conditions 
on the surface of the disk determined the scattering components. In 
order to insure that the tangential component of the electric field is 
zero on the edge, two sets of vector wave functions were used and their 
relative weights in the expansion of the incident plane wave were deter¬ 
mined by the edge boundary condition. This led to a unique solution. 
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The purpose of this work is to compare Meixner's and Flammer's 
solutions. Apparently, no numerical results have ever been obtained 
from Flammer's formal solution. Thus, it is necessary to establish that 
the formal solution is, in fact, valid. If this validity is established, 
the vector wave function formalism may be more convenient for examina¬ 
tion of such characteristics as the surface current distributions and 
the natural resonances. The reference work used for Meixner's solution 
is Hodge's version. Hodge[8] used the notation and normalization of the 
spheroidal functions introduced by Flammer. In order to match Hodge's 
geometry, Flammer's solution is rederived for the same configuration. 

The bistatic normal incidence case is more deeply studied, and problems 
encountered in the derivation of the solution are pointed out. Compu¬ 
tations made in the normal incidence case, as studied by Flammer, are 
compared to the results obtained by Hodge[9]. In addition, Meixner's 
fields are expressed in terms of spheroidal vector wave functions as 
well as the vector potentials. Another proof of Meixner's solution with 
vector wave functions is given. 
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CHAPTER II 

SPECIFICATIONS OF THE PROBLEM 


A. Geometry 

The scattering object is an infinitely thin, perfectly conducting 
circular disk of radius a. It lies in the x-y plane of a right-handed 
Cartesian coordinate system and is centered at the origin. 

The following study will be made in the oblate spheroidal coordi¬ 
nate system (n,£>4>), Figure 2.1, in which the disk is the surface repre¬ 
sented by £- 0 . 





1 < 0 


Figure 2.1. The oblate coordinate system. 


The transformation from oblate spheroidal to Cartesian coordinates 
are, as given by Flammer[10], 
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X 


'j-2a 


= a[(l-n 2 )0+f 2 )] 1/2 cos* 
jv - a[(l-n 2 )(lH 2 )] 1/2 sin* (1) 

[2 = a nC 


In Appendix A, the transformations between the different coordi¬ 
nate systems - Cartesian, cylindrical and spherical - and the oblate 
spheroidal coordinate system are summarized. Their limiting forms are 
listed for three cases: at large distances from the disk -£ large - on 
the surface of the disk and in the vicinity of the edge. 


B. Problem 


Because of the rotational symmetry about the z-axis it is possible, 
without loss of generality, to choose an arbitrary plane of incidence. 

We will use the notation of Hodge's work[9]. The incident plane wave 
direction lies in the x-z plane, coming from the x>0 half plane, Figure 
2.2. flammer used the y-z plane as plane of incidence, instead. The 
whole system is located in free-space. 



Figure 2.2. Geometry. 
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We will use the following symbols: 
c=ka - electrical circumference 
e 0 - angle of incidence 
a - polarization angle 
e 1u)t - time dependence 
The incident wave vector F is expressed by: 


F 1 = k[-sine 0 e x - cose 0 e z ] (2) 


The polarization of the incident electric field is measured b> 
the angle a between the plane of incidence and the incident electric 
field, Ei, in such a way that a is the conventional azimuthal angle, 
of a coordinate system defined by the unit vectors (cose 0 e x -sine 0 e z , 
iy, sine 0 e x +cose 0 e z ) as shown in Figure 2.3. Parallel polarization 
is obtained for a=0 and perpendicular polarization for a=-n/2. The 
incident electric field is given by 


tT 1 = E o (cos0 o cosa e x + sina iy - sin9 0 cosa e z ) 

• e' 1 ^ 1 ' r e 1 ' wt 


(3) 


where r is the usual position vector. 

Note that Flammer[7] characterizes the incident wave by the angle 
c between the positive direction of the propagation vector and the posi¬ 
tive z axis instead of the angle of incidence e 0 . Those angles are re¬ 
lated as follows: c=rr-e 0 . 

In the spheroidal coordinate systems, the boundary conditions sat¬ 
isfied by the electromagnetic field on the surface of the disk and at 
the edge are on the disk, z=0: 

(F + P) • e n = 0 

(F’+P) • e^ = 0 (4) 

(H’+H 5 ) • e r = 0 

where II s » E 3 is the scattered electromagnetic field. 
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Figure 2.3. Polarization 
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At the edge, c=0 and a=0, the energy density of the total field 
must be integrable. The components of the fields may become infinite 
at most as s' 1 / 2 , where s is the distance from the edge to the observa¬ 
tion point. Bouwhamp[10] showed the following more restrictive condi¬ 
tions on the ^-component of the fields: 


(r+£ s ) • e^ = 0(s 1 / 2 ) 

(HW) • e + = 0(1) (5) 

Those orders of variation as function of s are valid only in the 
vicinity of the edge, s<<a, for the curved edge of a disk. 

In this work, we will refer to these behaviors of the components 
of the total electromagnetic field in the vicinity of the edge as the 
edge condition. The ^-components will also be refered to as tangential 
to the edge, or parallel to the rim. 

We will see that, in Hammer's solution, the condition on [E 1 +E r ] ({j 
is equivalent to the whole edge condition. This is due to the particular 
choice of the vector wave functions used. This will not be true for 
Meixner's solution, where all the conditions on the components of the 
field will be needed. 

Before we present Hammer's and Meixner's solution we need to 
introduce the functions that we are going to work with and some of their 
properties. In Chapter III we will define the spheroidal functions, and 
in Chapter IV we will derive various plane wave expansions. In the fol¬ 
lowing work, the time dependence e iwt will be omitted. 
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CHAPTER III 

THE SPHEROIDAL FUNCTIONS 


In this chapter the scalar wave function solution of the scalar 
wave equation will be introduced and some properties of the spheroidal 
radial and angular functions will be summarized. The large radial argu¬ 
ment, c>\, approximation of the vector wave function solutions of the 
vector wave equation will be listed. 


A. Di fferential Equation and 

Scalar Wave Function 


In the oblate spheroidal coordinate system, the scalar wave equa¬ 
tion can be solved by the method of separation of variables. In Appendix 
B the different differential operators have been expanded in spheroidal 
coordinates. The scalar wave equation is: 


V 2 .j,+k 2 * = 0 

1 


2n + 1 -n 2 sfi + _2£_ + 1+£ 2 d P '<p 

r 2 +n 2 3n c 2 +n 2 8n 2 C 2 +n 2 n 2 +r 9£ 2 


4l- 


k 2 i it = 0 


The eigenfunctions associated with the eigenvalue x mn are: 


( 6 ) 


i / 1 ( i ) (n, £, <t>) 

mn 

o 


mn 


(-ic,n) 


r( 1 ) 

mn 


(-ic.ic) 



where S mn (-ic,n) is the angular function 


S«„Hc.r,> - y‘ d m "(-iO P^ (n) 


mn 


r=0,l 


m+r v 


as defined by FlammerpO]. 


(7) 


( 8 ) 
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P m (n) is the Legendre associated function as used in spherical 
coordinates. The expansion coefficients d l J! n (-ic) are computed from a 
recursion relationship as seen in Flammer's book[10]. 


The prime over the summation is a notation that will be used 
throughout the following work to denote that the summation must involve 
only the even or odd indexes when n-m is even or odd, respectively. 

- R^^(-ic,ic) is the radial function 

From Flammer's book[10]> 






2\m 


5" I <f n (-1c) &?^ r= °-’ 
r=0,l r 


l i r+fTW1 dj n (-ic) 


(2m+r) 
r! 


t-z ( i } (cO (9) 

n+r 


The index (i) denotes the kind of spherical Bessel function, z^ 1 ). 
used: i = l for j n , i=2 for n n , i=3 for h^ > and i=4 for h^). We will 
refer to i as the index of the radial function in the rest of this work. 
The normalization of the radial and angular spheroidal functions are 
those used by Flammer. 

In the following two sections, we will summarize some properties 
of the angular and radial spheroidal functions. 


B. The Angular Function S mn (-ic,n) 


The equation satisfied by the angular functions is: 


d_ 

dn 




-1 < n < 1 


The angular functions are orthogonal with weight function 1: 

rl 


-1 


^mn( _ ’ c »o)SmD(-ic,n)dn ^no ^mn 


J mp 


np 


( 10 ) 


(ID 


where 
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(cTHc )) 2 


n = ? V {r+2m)! _ 1 _ 

mn r=0 1 r! 2r+2m+l 


1. Evaluation at n=0 

From Formula (8), or more easily from the expansion of S mn (-ic,o) 
in powers of (1 -n 2 )» - Formula (3-2-7) in Reference [10] -, we can see 
that 

for n-m even S mil (-ic,n) is an even function of n 

for n-m odd S mn (-ic,n) is an odd function of n 


From the differential equation and its derivatives, we can evalu¬ 
ate the derivatives of S mn (-ic,n) for zero argument: 


2 [S mn (-ic,0)j “ ('\nn -m ) S mn (-ic*0) 


3 (S mn (-ic,0)) - (\nn~ m ~2) S mn (-ic,0) 


^~7T (5mn( _ i c ’0) ) -• [(^rnn _m2 ^ x mn“ m ^"^^ + 2m ~ 2c ^ 


Thus, the evaluation of S mn (-ic,n) and its derivatives around n=0 can 
then be described by taking into account the odd and even properties: 
n SMALL n«l 


W-Ic.n) = W‘ ic ’°) [l ‘ + 


S; n (-ic,n) = W-ic.O) [- nU mn -m 2 )] + 0(n 3 ) 


S mn(~ 1c > n ) 


- Smn(“ic,0) C - ( ^mn _rr|2 ) + 7 ~~ (^mn _rr * 2_ ^ 


+ 2m 2 -2c 2 ]] + 0(n 4 ) 
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where the normalization has defined S mn (-ic,0) = P™(0). 
n-m odd 

/ 

Snm(-i c »n) = S^, n (-ic,0)[n] + 0 (n 3 ) 

<S; n (-ic >n ) « S; n (-ic, 0 ) [1 - f- (x mn -m 2 -2)] + 0(n 4 ) (14) 

s mn^" ic » n ^ = s mn(~i c *°) C-n(A mn -m 2 -2)] + 0(n 3 ) 
where the normalization has defined S^ n (-ic,0) = ^[P™(0)] 


2. Evaluation at n=±1 

The power series expansion is led by (l-n 2 ) m//2 due to the P{{J +r (n) 
functions. For m greater than zero, 


m i 1 S mn (-ic,±l) = 0 
m > 3 S^ n (-ic,±1) = 0 


(15) 


The evaluation for m=0 and 1 has been given by Flammer[10]. 


S 0 „(-ic,l) * V d? n t-<c) 

r= 0 ,l 

s O n (-ic,l) = ~ l' r(r+l) d° n (-ic) 

L r=Q,l 

l im (l_n2) 1/2 i-[S ln (1c,r,)] -'-i I' (r+1)(r+2) dj. n (-ie) 
n+1 dr > m c r=0,l r 

SiJ-ic.n) ® In 

lim 21/2 * d r l*r+]0) 

n-*-l (1-n ) r=0,l 

= \ l (r+1) (r+2)d^, n (-ic) 

2 r=0,l r 


(16) 


C. The Radial Function Rj(^)(-ic,iO 


The equation satisfied by the radial functions is: 

2 


sr[( |w - 2 ) at - 1 


X r 2 f 2 
A mn c ^ 


in 


R^Hc.U) = 0, 


o < e 


(17) 


Unlike their spherical counterparts, Bessel, Neuman and Hankel 
functions - the oblate spheroidal radial functions and their derivatives 
are well defined over the whole range of values of f;. The normalization 
coefficients of these functions have been defined by Flammer[10] so that 
they match convenient formulas for large values of the argument c£, 



(18) 


With those definitions, we see that the suitable functions for 
the series expansion of an incoming plane wave will involve R^), and 
outward traveling waves will involve R^'. 


1. Wronskian 

The Wronskian of the differential equation is, as computed by 
Flammer[10]: 


a-.u, 


R^Hc.U) 


dR^(-ic,U) 
mn _ 

d£ 


= + 


1 

c(l+c 2 ) 


09) 


This can be written in terms of the R^, and R^ functions, using 
the property ^^R^-iR^), and the fact that the Wronskian of two 
linearly dependent solutions is zero: 
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t 




tfHc.U) ^ - R^Vic.ic) [R^Hc.ic)] 


c(1+c 2 ) 


(20) 


2. Evaluation at j=0 

In his book, Flammer lists the values of the radial functions and 
their derivatives at £=0. Appendix C of this work contains some correc¬ 
tions to those formulas. After presenting some general properties, we 
focus on R^) and RW separately. 

From the differential equation and its derivatives we can evaluate 
the derivatives of R^|\)(-ic,U) at £=0: 




CRjJhHc.iO) - (W<" 2 -2) 


Sun 


d£ 

& 


t R r in ,] (" ic » i0 ) = C(^mn” m2 )( x mn- m2 - 6 ) + 2m 2 -2c 2 ] f^(-ic,10) 

where c=ka (21) 


a. Properties of R^ 

Flammer[10] develops a power series expansion of R^P(-ic,i£). 
This function, if extended over the domain - 00 <£< 00 > would lead to an 
even or odd function for n-m even or odd, respectively. We can also 
conclude that, for k>0: 

7T = 0 (22) 

for k even or odd when n-m is odd or even, respectively. 
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In the vicinity of the disk, the radial functions of the first kind 
can be expanded as follows: 

f. Small c^ << 1 
n-m even 

/ 

= ^J ) (-ic,iO)[l+(A mn -m 2 ) |ij + 0(c 4 ) 

4 R mi ),( ' ic ’ i0 = ^J ) (‘ic.10)[€(A nm -m 2 )] + 0(c 3 ) (23) 

R il ),( "i C ’ ir ’) = R mJ )( ' ic » i0) | A mn- m2 + C(W" 2 - 6 > + 2m 2 -2c 2 )] 

^ L + or^i J 


n-m odd 


^^(-ic.iO = ^'(-ic.iO)^] + 0(c 3 ) 

R^’Hc.U) = Rin r (-ic,i0) ^l+( A mn -m 2 -2)§^ + 0(C 4 ) 
(^^"(-ic.ic) = Rj 1 1 n ), (-ic,iO).-(A mn -m 2 -2) + 0(c 3 ) 


b. Properties of 

From Equation (19), the Wronskian of rI^ and R^ 4 ^ leads us to the 
following identities 


n-m even R< 4) ' He,TO) = - --(V . 


c C (- ic >-°) 


n-m odd R^(-ic.iO) = -pryr 1 - 

c ii ] n } C-ic.iO) 


We then get the following expansions in the vicinity of the disk: 
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For c Small, cf << 1 


n-m even 


R^i n ^(-i c >ic) ~ l4 n ^(-ic,i0) |l + U mn -m 2 )j 


M7~. ~ + 

c ^mn ( - ic,iO) 


+ ^ (A nrn- m2 > f 4n ) <- ic » i °) + OU 3 ) 

^"Hc.U) = ^^(-ic.iO) (26) 

[ (A mn- m2 ) + ^ A mn-m 2 )(^ n -m 2 -6) + 2m2-2c 2 )]J 

u c R^-icJO) C(x mn' m2 - 2 ) + OU 3 ) 

Note that the coefficients of R^-ic.iO) in these equations are 
( 23 ). t0 th ° Se ° f ^n (~ ic »i0) in the corresponding equations of Formula 


n-m odd 


R^Mc.U) - —7-| \ - i - 1 + ( A mn‘ m2 )] 

c (-ic,iO) L 2 m " J 

+ 5 Rm^'Hc.io) + OU 3 ) 

R mS ), (- ic ’ i ^) = R in > ’(- 1c - , 0) [ + I" ( A mn' m2 - 2 )] 


+ ~jrrrr~7^ ^w* 2 ) + °u 3 ) 

c Km Hc.iO) 


Hc,u> ‘ r^jV(- ic .(o) [ in, "'" z + r- 

♦ 2B 2 -2c 2 ]j + 5 R<J) , (.i ti io)(» |nn V-2) * 0(t 3 ) 
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Note that the coefficients of (-ic,io) in these equations 

are equal to those of (-ic,iO) in the corresponding equations of 
Formula (24). 

3. Evaluation of R^ for large arguments 
Equation (18) gives the large argument approximations 


knn*H c »U) 


exp(-icf) 




¥ [■• • i] 


exp(-ic^) 


[R^Jht-ic.ie) - > [= 2 - ~ - ^ ex PHcc) 

b 

In the far field expansions, only the terms with dependence 
will be kept. 


D. The Vector Wave Functions 

Solenoidal solutions of the vector wave equation can be expanded 
as summations of spheroidal vector wave functions. Those will be de¬ 
fined, following Flammer's notation, by: 


M e “' 1 '(n»f ,<t>) 


= v x [^ i; (n,C»*)e a ] 


N p a ( i )(n,f,4.) 

omn 


I v x M! (i) (n,C,<!>) 

K v mn 


where e is either a constant vector in Cartesian coordinates (a=x,y,z) 
or the position vector f=r S r (a=r). 

Extensive work has been done on these functions by Flammer, who 
lists them all in his book[10]. In contrast with the orthogonality of 
the scalar wave functions, the spheroidal vector wave functions are not 
orthogonal. They are not even independent as will be seen in Chapter 
IV-C. We will list here the behavior of those functions on the surface 
of the disk and for large radial arguments. 





1. Behavior of the vector wave 
functions at f,=Q 


We have seen Chapter III-C, Equation (21), that, on the surface 
of the disk, all radial functions and their derivatives can be expressed 

in terms of (-ic,iO) and R^i) (-ic,i0). We can show that, on the 
surface of the disk, each component of any vector wave function defined 
by Equation (29) is proportional to only one of these values. We list 
the different cases. 


a. Tangential components at g=0 


[M x(i)] , [mX( 0 ] , [M* (i) ] , [M* (i) ] ,[F£ (i) ] and 

<ynn n e mn < p e mn i e mn <P e mn n 

v 0 0 0 0 

are 
r>mn 


/ i \ I 

proportional to R^ n ' (-ic,iO). 


(30) 


:;mn 


• [ SX e‘ i> ]* • • [ ^’ ] n • 

^mn ^mn * ^mn Qmn v ^ In 1 

are proportional to R^(-ic,iO). (31) 


b. Normal components at C=0 


[N e x(i) L 
e mn 5 

0 

, [N y(i)] 
o mn 

are proportional to R^ (-ic,iO) 

(32) 


> [M/ (i >], 

[N e z ^] are proportional to 


o mn ? 

0 mn £ - 

„mn £ , 

(33) 


We notice that at the surface of the disk the tangential components 
of any of these vector wave functions, and_the normal component of their 
curl - that is, for example, [N 01 ]^ ^ and [M a ] £ or vice-versa - have the 
same radial dependence. That ensures that both boundary conditions on 
E and H at the surface of the disk are matched by the same expansion. 
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2. Behavior of the vector wave functions 
of index i=4 for large arguments of 
the radial part 

In this section, we list the large radial argument approximations 
of the vector wave functions, keeping only the 1/f terms. From Equation 
(29), 


f M e X(4) = f f- S m n ( - ic >n ){ - sin 4> e n +n cos* e^} ™*m* e' 1cC 
0 mn 


V (4) = | 7“ S mn (-ic,n){cos<j» e n +n sin* e^} £0%* e -icc 

r>mn 


\ '^1-n 2 S mn (-ic,n) ^ r^m* e e H 


e m „ 

0 mn 


coSfpA «-i cz 
4 sin 


(34) 




,-n-l 


cos* e n +sin* e^} ^m* e~ lC ’ 


N|*‘ 4) ■ !r~ W-fc.n)U 
ff y(4) s 2 i/^l s mn (-ic,n)(n sin* e n -cos* e^> c ^ s m* e' lcC 


e H 7- mn 

n mn d f 


sin 


z(4) * _ <■ /f. n 2 S mn (-ic,n) cos m* e e~ ic S 

mn n 1,1 • * r- ~ A n 


-mn 


t, sin 


We now have enough tools to find the expansions of scalar and 
vector plane waves in terms of the spheroidal functions. 
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CHAPTER IV 

PLANE WAVE EXPANSION 


In order to expand the incident fields, we need the expansions 
of vector plane waves in terms of the spheroidal vector wave functions. 
First, we will derive the scalar plane wave expansion in terms of scalar 
wave functions. Taking its curl associated with different basis vectors 
will lead us to various vector plane wave expansions. The special cases 
of normal incidence to the disk are reviewed. 


A. Scalar Plane Wave 


By scalar plane wave we denote the function exp[ik(x sino Q + 
z cosGg)] in our geometry. We do not know enough about the angular 
spheroidal functions to directly obtain the coefficients of the expan¬ 
sion of the scalar plane wave in terms of spheroidal functions by 
applying orthogonality on the trigonometric and angular functions. We 
will use another method as derived by Flammer[10]. The plane wave is 
considered to be due to a point source removed to infinity in the 
direction of arrival of the incident wave. The scalar plane wave will 
then be obtained from the asymptotic form of the free-space Green's 
function. 

Because of our choice of the time dependence, e iu)t , the Green's 
function used here will be the complex conjugate of that obtained by 
Flammer[10]: 

exp(-ik|r'-r*I) 

4tt | r'-r' | 


(35) 


CO 03 


lk i_ I nr ? 21 


2tt 


m=0 n=m ron 


R^(-ic,U' )Rmn^(-ic»iC) cosm(<t>-<^') when £*>?, 


where (n',f',<f>') are the coordinates of the source point. 

For a point source at infinity, we can use the approximations: 


exp(-ik|r-r'I) m 
4tt| r-r‘ | c£'-*» 


exp(-ik(r'-r.y l r l i'| 

4<i r' 


(36) 
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R iiS ) (-ic,U‘) 


c r 


±r exp;-i(cc' - ^ -)J = ^ e" 1 ^ 


cc* *-•* kr 1 as seen in Appendix B 

c 1 -*» 


Finally 


exp(ikr*v'r') = 


2 Jo l ^T“ 4n ) (- ic .C0S9 0 )^’ ) (-fc,n) 

x 4^(-ic,i'J cos m U-*’) (39) 


r-V'r' - r *[si nO Q cos<t>' e x + sine 0 sin^' e y + cos9 0 e 2 ] 
fVr' = x sino 0 cos<{> 1 + y sine 0 sin4>‘ + z cose 0 


For a wave incident from the positive x half x-z plane, <»'=0 
and we obtain the expansion 

* 00 2-On m 

exp[ik(x sine 0+z coso 0 )] = 2 [ I lj—- i S (-ic.cose ) 

m=0 n=0 mn 


Smn^ c,n )*4n^'i c »"'^ cos m$ . (41) 


We define the coefficients, Y mn (° 0 )» of the scalar wave functions, 
as in Flammer[10], Equation (7): 

(2 -‘W . 

' r mn'°o^ ~ ^ N mn ^ ^mn^“^ c,coso o^ (42) 


•r> r,t / 

exp[ik(x sino Q + z coso 0 )] = \ \ Y mn (8 )*^ ] (n,C ,*) 

m=0 n=m mn 


2 Q 



B - Vector Plane Wave Expansion 

■“■ara.'isaa svr- 

(29). ^ liSt thG obtained formu las, with the notations of Equation 


e y exp[ik(x sine Q + z cose 0 )] 


L 1 Y mn (f U M P x(1) -n,f,o) 


ik cose o A n=m W'V \ 


(-e x cose 0 + e 2 sine 0 )exp[ik(x sine 0 + z cose 0 )] 

-j 00 oo 

Tk ^ Y mn^ e o^ ^emi 1 

m=0 n=m mn 

e y exp[ik(x sine Q + z cose 0 )] 


ik sine o m=0 n=m ){n,C><(>) 


( _e x C0S9 0 + e z sine 0 )exp[ik(x sine 0 + z cose 0 )] 

<*> <x> 

5 - r~r L 1 W 8 o) R e!!„ (■'■«■♦) M 

K COS U Q m= Q n=m C mn 

A OO OO 

e J exptutx S ine 0 ♦ z cose 0 )] - l J l v„, n (e 0 ) il/' 1 

m=0 n=m e mn 

(4 

(" e x cos9 o + ®z Sin9 0 )exp[ik(x sine 0 + z cose Q )] 

" rik; J 0 l Y »" <e ° > N em 9)(4! 
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C. Normal Incidence Case, o 0 =0 

In this section, we will study the behavior of Equations (44) to 
(49) in the case of a plane wave normally incident on the top of the 
disk. As we have seen in Formula (15), the angular functions of order 
greater than zero are equal to zero for unity argument. 

m£l S mn (-ic,cose 0 ) = 0 for e 0 = 0. 

From Equation (42), we deduce that yg n (0) is the only non-zero 
expansion coefficient. The expansions Equations (44), (45), (47) and 
(48) can be readily seen to be single summations on n, 0<n<«> with m=0. 
The expansions Equations (46) and (49), however, involve a l/sine 0 
factor that leads to indeterminacy of the limit when e 0 tends towards 
zero. 

As we can see from the definition of S mn (-ic,n), Equation (8), 
the angular function is led by a (1- n 2)m/2 factor. Therefore, from 
Equation (42), we can write: 

YmnK) “ O-cos 2 e 0 ) m/2 = sin m e Q (50) 
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For m=0, YQ n (o 0 )/sine (;) goes to infinity when e 0 approaches zero. 

We will show that, despite the fact that each coefficient diverges, the 
partial summations over n for m=0 in Equations (46) and (49) are globally 
equal to zero. This will lead to two dependence relationships for the 
spheroidal vector wave functions. Consider first the case of Equation 
(46). Because of an m coefficient, the n- and ^-component of M e l^' are 
equal to zero. We prove that the ^-component of the partial summStion 
for m=0 in Equation (46) is also equal to zero; this is expressed in the 
following equation: 



TP S 0n (-ic,n) R^(-ic,U)-£ S 0n (-ic,n) (52) 


<Lr 
de K 0n 


{1) (-ic,ic)]=0 


'J 


This must be valid for any pair (n»0- In order to prove this 
identity, we must go back to the scalar plane wave, Equation (41), and 
specialize to our normal incidence case. 


exp[ikz] 


°° (11 

I W°) W' ic ’ n ) ^n Hc.ic) 


(53) 


From Appendix B, we know that z can be expressed very simply in 
terms of n and p 

z = an? . (54) 

Taking the derivative of Equation (53) with respect to n and r, 
separately, we obtain two expansions as follows 

oo 

~ [exp(ikz)] = ikac exp(ikz) = l Y 0n^ 0 ^ S 6n^“ ic,n ^ R 0n^" ic,i ^ 
an n=0 

(55) 


fr [exp(ikz)] = ika n exp(ikz) = l y Qn (0) S Qn (-ic,n) R^ 
^ n=0 


(-ic.ic) 

(56) 


Multiplying Equation (55) by n and Equation (56) by -c and adding 
them leads directly to Formula (52). We have shown that, for any Dair 
(n»6): 





5 


( 57 ) 
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Taking the curl of this equation leads to: 


Jo W 01 S 4T !( "- f -* ) ' 6 (58) 

To suppress the indeterminacy of the limit of the ratio of Equations 
(57) and (58) to sine Q when e 0 goes to zero, we will apply L'Hospital's 
rule to each component of those equations: 


lim -J— l y (e 0 )[FL 

.n sine„ nEn On ° /L 8 


o >0 sine 0 n=0 


= Hm [ 


?i n S 0n (ic,cose 0 )(-sine o ) 

- 1 —- ■ — 1 - --- r n 


o 0 >0 n=0 "On 


where [F] e represents the B-component of either M p z ^ or N z ^) 


8=n»t, 


T 5 U //,nt conc J ude that > in the normal incidence case, the expansions 
(46) and (49) can be expressed as single summations on n, l<n«= with m=l 
[;lammer[7] gives a vectorial proof of Equation (57) by taking the curl of 
e z exp(ikz). Equation (53), which is zero as the curl of a z-directed 
vector whose component is a function of z only. 

For convenience in the study of the normal incidence case, we de¬ 
fine, as FTammer[10], the following coefficients: 

oo 

-On = 2 *"■' < V d? n (-ic) 
r=0,l 

■ 2 K 0n 5 0n (-ic,l) 


I 

i Y 0n 




b, * 2 1" +1 n: t T' Lam d Tn (-i C ) 

In In r 4_, r! r < ’<=) 

The vector wave expansions are then, for normal incidence: 


e exp(ikz) = i f a M x(1 
y k n=Q On e Qn 
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( 63 ) 




In the previous two sections, we have derived various relation¬ 
ships between spheroidal scalar functions, and between spheroidal vector 
wave functions. Some limiting cases will allow us to show some depend¬ 
ence relationships, as follows. 


Equations (55) for £=0 shows that the derivatives of S Q (-ic,n) 
are not independent and therefore not orthogonal as the original func¬ 
tions were. 


Equations (57) and (58) show the lack of independence of the func¬ 
tions M 6 Q n , and N 6 q^ ^, respectively. Equations (44) and (47) ^multi¬ 
plied by cose 0 , lead to dependence relationships for M e x O), and N g x ( ' 

respectively, when 0 g approaches tt/ 2. That kind of rela?ionship can n be 
derived for each family of vector wave functions of index i=l. We do 
not deal with independent sets as in the spherical case. In the next 
two parts, we will derive two different solutions for the scattering 
problem. In both cases, the fields will be expanded in spheroidal vec¬ 
tor wave functions. For the incident fields, we will use the expansions 
derived in the previous two sections. We will not be able to consider 
individual terms in the summation expressing the incident field to be 
modes as we do in simpler problems, where orthogonality properties exist. 
To each vector wave function of the expansion of the incident field will 
correspond not only one vector wave function in the scattered field but 
a summation of the scattered counterparts of all the vector wave func¬ 
tions appearing in the dependence relationship satisfied by the initial 
wave function. This will be shown in part VI-C. We have not, however, 
proved that the spheroidal vector wave function of index i=4 are 
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dependent since the relationships obtained for index 1 are not appl 
cable to functions of index i=4. 

In the following chapter, we will rederive Flammer's solution 
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CHAPTER V 

FLAMMER'S SOLUTION 


Unlike Meixner, F1ammer[7] deals directly with the electric and 
magnetic fields in order to find a solution to the scattering problem. 
This approach allows simpler, more physically interpretable boundary 
conditions, but at the expense of more complicated expansions. The 
major problem encountered in the scattering by a disk is that satis¬ 
fying the boundary conditions on the surface of the disk leads to solu¬ 
tions that do not necessarily satisfy the edge condition. The edge con¬ 
dition puts restrictions on the behavior of the components of the total 
electromagnetic field in the vicinity of the edge as explained in Chap¬ 
ter II of this work. For the electromagnetic energy density to be inte¬ 
grate in the vicinity of the edge all components of the fields vary at 
most as s" ' where s is the distance from the edge, while the ^-compon¬ 
ent of the total magnetic field remains finite and the <(>-component of 
the total electric field goes to zero as s'' . A very remarkable point 
in Flammer's solution is that the satisfaction of the condition on the 
(fi-component of the electric field directly ensures the satisfaction of 
the edge condition. 

Whatever the expansion of the electric field of the incident plane 
wave in spheroidal vector wave functions of index 1, the summation of 
spheroidal vector wave functions of index 4 which satisfies the E-field 
boundary conditions on the surface of the disk has a singularity of 
order s'^ i n its ^-component. Flammer chose To express the incident 
plane wave as sum of two different expansions. The relative weights of 
each expansion are then chosen in such a way that the coefficients of 
the singularity in the ^-component of the corresponding E-field expan¬ 
sions of index 4 cancel each other. With the proper choice of the two 
plane wave expansions, the above procedure also leads to the satisfac¬ 
tion of the whole edge condition. 

First we will study the normal incidence, perpendicular polariza¬ 
tion bistatic case because of its greater simplicity. Perpendicular 
polarization is here understood as defined in Chapter II, considering 
the x-z plane as plane of incidence. The problems encountered in ex¬ 
tending the method to the arbitrary incidence, arbitrary polarization 
bistatic case will be noted, and the results of a numerical test of the 
far-field for the normal incidence case will be presented. 
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A. 


Normal Incidence, Perpendicular 
Polarization Bistatic Case 


This^ corresponds to the case a=ir/2, £P)|e y , and FP||e„. The 
incident E-field can be expanded according to Equations (62;, (64) or 
(66) and the corresponding H-field according to Equations (65), (67) 
or (63), respectively. To represent the outgoing scattered waves, we 
use vector wave functions of index 4 which match spherical waves when 
t. becomes large. 

We will successively consider the boundary conditions on the sur¬ 
face of the disk, an expansion of the ^-component of the electric field 
in the vicinity of the edge, the [P], condition, and the final solution 
where we will check the satisfaction^ the edge condition. 

1. F ield boundary conditions 
on the surface of the disk 


We will successively consider the three possible expansions and 
the scattered fields they lead to. Despite the lack of orthogonality 
of the spheroidal vector wave functions of index 1, we will use a term 
by term, function by function, matching technique to obtain the expan¬ 
sion of spheroidal vector wave functions of order 4 which satisfies 
the boundary condition on the surface of the disk for given incident 
waveexpansion. These expansions will be called "reflected" fields P" 
and H r . The reflected fields are not equal.to the scattered fields 
unless the edge condition is satisfied by E’’+P' and H 1 +H r , but they 
satisfy the boundary conditions on the surface of the disk. The bound 
ary condition is given in Equation (4): 

[f+E 5 ], = [f+E 5 ]^ = [hVh 5 ]*- = 0 


The term by term matching technique does not lead directly to the scat¬ 
tered field since the spheroidal vector wave functions of index 1 are 
not independent. We now express the reflected fields for the three 
possible expansions of the incident fields. 


- x( i ) 

M q v ' expansion 
_ 

From Equations (62) and (65), 
E 


FI . 1 V, m x O) 
' k n = 0 0n e 0n 


n 1 - i H ° y a N x ( 1 ) 
H - ' r „i 0 a 0n %„ 


where H 0 =E 0 /Z 0 and Z 0 is the free space impedance. 


( 68 ) 


i 

t 
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From Equations (30) to (33), and Equation (22), we see that, for 

n even, the components of M *0) tangential to the disk and the normal 
_ J(d e Qn 

component of N e ^ ' are equal to zero. Therefore, in Equation (68), 

the vector wave functions for n even do not lead to any reflected.field. 
For n odd, the above three components are all proportional to RX^(-ic, 
iO), and the same components of the functions of index 4 are propor¬ 
tional to R^J)(-ic,iO). Since the only difference between a spheroidal 
vector wave function of index 1 and its counterpart of index 4 is their 
radial dependence, it is possible to cancel a component of one by that 
of the other on the surface of the disk with appropriate coefficients 
in order to satisfy the boundary condition. 


F( ] )(n,CM) 


R 0n ( - ic > i0) 
Rg^ (- i c, i 0) 


F^ 4 ) (n,0,<!>) = 0 for 0<|n|lO 


where 


F^ )( n £ <t>) = I'M X ^1 FM h FN x(i) l 

i- Li\ 0n J n » L^e 0n \ » L"e 0n h 


The term by term cancellation method leads to the following 
reflected fields which satisfy the boundary conditions on the surface 
of the disk 


ft _ Si y . R 0 n 1 ' 1 -’ c - ,0) jfX( 4 ) 

k n-1 °" R‘J>'(-ic,iO) e °" 


RIP (-ic,iO) _ 


nr = _o yi On __ 

k rJl 0n R^'f-iciO) e °" 


where the prime over the summation means that the index of the summation 
varies by increments of 2 from its initial value. This notation will be 
used throughout the rest of this work. In the same way, for the other 
two expansions, we obtain the following formulations. 


y*(D 


expansion 


From Equations (64) and (67), 
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( 72 ) 




From Equations (30) to (33), and Equation (22), we see that for n-1 
odd, the components of ^ tangential to the disk and the normal 

component of are equal to zero. For n-1 even, the same com¬ 

ponents are proportional to Rj^(-ic,i0) and the corresponding compon¬ 
ents of the functions of index 4 to r(^(- ic,iO). 

I n 

Using the same term by term matching technique, we obtain the 
following reflected fields: 

E. Rp^(-ic,iO) _ 7 ( 4 ) 

ft Ori, In m z ' ' (7 n' 

E - ■ r b] " ( ' 


i H '» 

TT I; b ln 
n=l 


Rfl^-ic.iO) _ z( 4 ) 

—rr\ - N 


R)n ("iC’iO) 


- y(l) 

No expansion 
K 0 n _ 

From Equations ( 66 ) and (63), 


^ "4- I a 0„ . ^-rl a 0n * 


m7(D 


k n^O °" 6 0n 


^ n=0 d ° n e 0n 


From Equations (30) and (33), and Equation (22), we can see that for 

n odd, the components of ^ tangential to the disk and the normal 
_ v (i ) e 0 n 

component of M--; are equal to zero. For n even, the same compon- 
e On 

ents are proportional to Rg^ ‘* 0 , and the corresoonding components 

of the functions of index 4 a>- e proportional to Rg^ 4 '(-ic,i0). 

Using the term by term matching technique, we obtain: 


E = - 


P a R 0 n )( ~' C - i °) s y(4) 
k n=0 0n R^)(-ic,i0) e ° n 


30 


H 


(77) 


k n=0 0n R^I-ic.iO) 


m y(4) 

e 0n 


All three pairs of incident and reflected fields satisfy the 
boundary condition on the disk, excluding the edge. We will study the 
behavior of [E 1 +F r ] j j ) at the edge in the following section for each case. 


2. [E]^ condition at the edge 

We will first give the expansion of [E 1 +F r ] ( . in the vicinity of 
the edge, and then present Flammer's method. 


a. Behavior of [F+eH 

___i 

at the edge 

The edge is located at the coordinates n=0 and £=0. We will use 
the small argument expansions for the.angular and radial functions to 
obtain a power series expansion of [E i +E r ] ( j ) . In each summation, we 
will only retain the diverging and constant terms, when they exist, and 
the converging term of smallest order at the edge. We will use the 
following expansions: 


(l-n 2 ) 1/2 =1 . 1 n 2 + 0(n 4 ) for n and £ small 

d- n 2)-l/2 = t + ^ n2 + 0(n 4) J n |« 1 

( 1 h 2)!/2 = 1 + ^ c 2 + 0 (c 4 ) e « i 

( 1+£ 2 )' 1/2 = 1 - 1 £ 2 + 0 (£ 4 ) 


We will now consider the three cases. 


— x(l) 

Case of the M e ^ expansion 


(78) 


A 


For this case only we will give the whole procedure: 


t 
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[ £l+E, % ; k )•' a 0n d 2 2 (1 - n 2 ) n (- a Qn ) ,1 + Ag n j -^ 

n=0 +n L \ 


♦ 


2 “ On / 


. I y a 2 _J_ 

k i a °n d 2 2 

K n=l u f 3 +n 


+ n(l+C 2 )n 6 + §- ( ; 0n 


jl v 1 , - un 

k Z a On d -T4 
n=l 


n(l+C 2 k >- 0n j Ro ^( "ic>iO)S 0n (-ic,0)cos^ 
-J -2 [f.O-n 2 )(l U 0 n- 2 >)£ 

+ §“ * ; On' 2 y- R &n*'H c ’ i0 ) s 0r/~ ic '°) coss 

' 

h 2) KJ 


1 V, 2 R^McJO) 1 f n 2 


R On Hc.iO) ^ 2+r|2 


‘ 2 ^ A On~2)) 


+ nO +r /)n£ A, 


C(1V) 


i S' (-ic,0) 


c RQ n ' (-ic.iO) 


+ C(l-n 2 ) 


(l - r f~ ( A 0 n _2 ^ 


C+n(l+C 2 ) n(l +|-(>fl n -2) 


RqJ } (-ic,i0) x S' (-ic,0) 1 cos? 


We can see directly that the n=odd terms of the incident wave are 
cancelled by part of the terms of the "reflected" field in this expan¬ 
sion. Finally Equation (79) reduces to: 

[E 1 +E r ] * jj ns V a 0n » 0n ^^’(-‘c.lOJSonf-tc.O) 

n—u i 


n 2+ ^, 2 n =l c 


f, t (n 2 n 2)] 

Ri 4 


R 0n >,( - iC * i0) | 

X COScf) (80) 


The singularity at the edge comes from the second term, 
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r r, 1 "'On 
„ V n=l C a ° n ^ 


Snn(-ic.O) 


R On ( _i c,iO) 


which does not converge when and n->-0. 

Case of the Mp^ 1 ) expansion 
_Ijj]___ 

The same procedure, taking into account the use of the approxima¬ 
tions of Equation (78) leads to the following result. 


■ - h 


A ln' 3 




n € l b ln-l- S i n ("ic*0)Rln Hc.iO) 

n=2 


+ -f-o X'i -la i + J" (n 2 +c 2 } in - 

nV n=l c L 2 J R|n ; (-ic,i 


0) 


COSO 


(81) 


We notice that the singularities in Equations (83) and (81) are 
of the same order C/(n 2 +^2). This is the reason why Flammer is able to 
cancel the singularity by taking the incident field as a linear combi¬ 
nation of Equations (68) and (72) with proper relative weights of the 
two expansions. 


Case of the 

___^On 


expansion 


The [N e J (1) ] 

angular and raOial 
compute. We will 


^ component involves second order derivatives of the 

functions and are, therefore, more complicated to 
use the following expression for that component: 


[N y(0] 

e 0n « 


* ' W 2 jf? ^ (s ° n> 

+ (, - n2 > <W ^ < R 0n’> 


+ <'* ?2 > s 0n 



(82) 


Using Equations (13), (14), (23), (24), (26) and (27), we finally 
obtain for the tangential electric field at the rim: 
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[F+r'i 


k 2 d 2 


cost 


■: 3(* 0n -2) 

+ lrr ' Ll a 0n -2- s 0n(- ic >°) 

n= I 


m ' 

Ron Hc,iO) + 


n 2 +f 2 n=0 - 


I 1 + * 2 (-o„* 2 > - 


(A 0n -2) 


a Qn W" 1 ^ 0 ) 


c Rfln ^ (' 1 c > i 0) 


( 83 ) 


These three expansions of the ^-component of E’i+E r for different 
incident plane wave expansions behave in the same way in the vicinity 
of the edge. They each have a singularity which comes from the ill-de¬ 
fined spheroidal vector wave functions at the edge. We have seen by 
comparing the power series expansions of r( j/(-ic,ir) and its derivatives 

for small argument Equations (23) and (24), to those of R^'(-ic,ic), 
Equations (26) and (27), that part of the expansions of the TOnctions of 
index 4 is identical to the whole expansion of the function of index 1 
except for the change of index. The term by term matching technique 
exactly cancels these terms between the incident and reflected fields. 
Therefore, in Equations (80), (81) and (83), the second summation con¬ 
taining the singularity comes from the terms in the expansions of the 
radial functions of index 4 that do not disappear in the term bv term 
matching. Their coefficient r/(o 2 +f 2 ) shows that, as expected, EC i +E r is 
equal to zero on the surface of the disk. Furthermore, the first sum¬ 
mation in Equations (80), (81) and (83) corresponds to the part of the 
incident plane wave expansion that does not lead to a reflected field 
since it is equal to zero on the surface of the disk. This summation is 
also equal to zero at the rim as the factor ru; shows. 


b. Flammer's method 

In order to satisfy the [E]. condition, Flammer[7] simply chose 
to express the incident wave as the sum of two different expansions 
whose weights would be adjusted so that the singularity in the compon¬ 
ent of the electric field at the edge would disappear. In order to 
satisfy the whole edge condition by the procedure, it is necessary to 
choose the two expansions so that the other components of the fields 
haye the right behavior near the edge. Tie n- and t-components of 

^ have higher order singularities than those of M p *' ' and 
e 0n 3 3 e 0n e ln 

while the latter have the same singularities, and so do their magnetic 
field counterparts N e ^ ] ) and N^ 1 ). We will show that the choice of 
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the two expansions in terms of M Equation (68), and 

e 0n e ln 

Equation (72), leads to a satisfactory solution. 


1 a °n< 1} +3 2 I hn^ 1} 


with 3 X +B Z =1 in order to have |F|=E 0 . 

The scattered field will be: 

n 1 • 

^ E o „ R 0n Mc.10) si x (4) 

E " k P 7 a °" R^'Hcio) e °" 

°° R ln^("ic»iO) _ 

+ l' b ln -7 M e? (4) ' (85) 

n=l R ln (“iC’iO) ln 

From Equations (80) and (81), we see that the coefficient of the 
singularity of the expansion of [F+E s ]^ is: 

o “ -i s n n ( _ic »°) „ V' 1. s i n (' 1c ’°) 

— _ o Y 1 — a* i i . -+ 3z £ p b-i m J C0S( p 

kd B x t c a 0n 7(4) . .77 n=l c ln F 4 


R 0n ( - ic * i0 ) 


(-ic,i0) 


In order to set this coefficient equal to zero, we must choose 
3 X and 3 Z as follows: 



With b x and 6 defined as above, F+F satisfies the [EF condi¬ 
tion at the edge. Equation (85) gives, then, the bistatic scattered 
field of.the disk with normal incidence and perpendicular £olarization 
(i.e., F||e y ). In the next section, we will verify that E s satisfies 
the whole edge condition, underlining the convenience of Flammer's met¬ 
hod which satisfies the edge condition by dealing only with the condi¬ 
tion on the <i>-component of the electric field. 
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3. The bistatic scattered field 


The total scattered electromagnetic field obtained with Flammer's 
solution is, in the normal incidence case: 


= 


H 5 


\< i: 


n=l 


(’ll 1 

R 0n Hc.iO) - x(4) 

a °" i^HcVio) 


x - Rp^-ic, iO) 

+ [i z l' b ln ~m-~~ M ei 

n=1 R ln (” ic *10) ln 


iH r 


< i :* 

n=l 


fir 

R 0n (~i c » i0 ) — x(4) 


On JJT 


Rgn' Hc.iO) 


e 0n 


+ y. b R 1n >( -' ic - 1 ° > „ zO> 

z n=l ln R[^(-ic,iO) 6ln 


where and are defined by Equation (87). 


( 88 ) 


First we will verify that the above fields satisfy the whole edge 
condition in the vicinity of the edge. The behavior of each component 
of the fields will be expressed in terms of local coordinates that have 
been defined in Appendix A. Second, the bistatic scattered far field 
will be computed in spherical coordinates. 


a. Total field in the vicinity of the edge 

The edge condition is expressed in terms of the distance from the 
edge to the observation point. We will introduce the local coordinate 
system, valid for n and T small, shown in Figure 5-1. We use the follow¬ 
ing notation: 

s = distance from the edge, s<<a 

t = angle between the top surface of the disk, n>0, 
and the direction from the edge to the observation 
point, 0<t<2n. 

I = same as in spheroidal coordinates. 
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Figure 5-1. (s,t,<j>) coordinate system. 


The 

ordinates 


correspondance between spheroidal coordinates and these co 
is derived in Appendix A. 


(s = j (n 2 +? 2 ) 

4 


2 .2 

cost = - 2-2 
n +e 


n +C 


or, conversely. 


(89) 




sin(t/2) 


(90) 


The disk is characterized by t=0. 

In order to obtain the components of the fields in the new coordi 
nate system, we first expand the radial and angular functions using 
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their small argument approximations and then transform the terms of lower 
order with the above formulas. 


For the n-component of E=E 1 +E r , the singularity in [M e *^] n and 

[M e ^^] n is (n 2+ T 2 )"^ 2 , or s - ^/ 2 , and [lf] n has a satisfactory behavior 
near the edge. More precisely, we have 


= E o ej r l' a 0n (-!) sint/2 S 0n(-lc,0)R<n , '(-ic.i0)s) 
* n=0 x f 

“ (1) 1 

a ° n (-j)^ s1nts »" ( - ,c,0> -zrS 


sint S on^' ic, °^ 


RqJ^ (ic.iO) 
R 0n^' Hc.iO) 


x r( 4) (-ic,iO)sind) 

On 

+ E 0 3 z l J’ b ln I JI sint S ln (-ic,0) 

x (-ic.iO)sin<}> (91) 

In 

The first summation is equal to zero on the surface of the disk 
but apparently only there. However, if we differentiate Equation (56) 
again by £;, we obtain 

2 °° 

—o [exp(ikz)]= -k 2 a 2 n 2 exp(ikz) = l 7 nn (0) S Q (-ic,r.) 
ar/ n=0 un 

' 1)" 

x R 0 n Hc.U) 

(92) 

At the edge, where f,=0 and n=0, Equation (92) becomes: 

).'• T 0n (°) S 0n (-ic,0)R^ ) "(-ic,i0) = 0 . (93) 

n=0 


R](~ic,i0) 
R]|P(-ic,i0) 
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The terms of this summation for n odd disappear since, for n odd, 
SQ n (-ic,0) is equal to zero, Equation (14). Equation (93) is the coef¬ 
ficient of sint/2 in the first summation of Equation (91), recalling 
that by definition a 0 n =-iAQ n (O)• Finally we see that Equation (91) 

contains no term larger than s ^/2 and therefore, we have at the edge: 
[E] n = 0(5^2) sin<f> sint (94) 

For the e-component of E, we see that the singularity of 

and [M -j s due to a f ac tor (l/(n2+e^)^)» or s"^2 > u ee the 

definitiBn of the functions in Flammer[10]) and therefore [FL satisfy 
the edge condition for this component. 


[E] f = E 0 I' a 0n f cost/2 S 0 n (-ic,0) Rj n ; (-ic,iO)sin* 


01 “ n a nn SA n (-1c.O) t 

+ E x y. l_0n On - LJl sin* 


+ E X y. 1 _yn - IdLL sin* 

° k n=l c ^ rJJ' (-ic.iO) /s 

8 _ r 4i r~ 7 ^"ln(“i c >^) 

+ E o ir ^ b in%72 /? sinm ~vrr, -r sin * («) 

k n=l cd J/ 2 Rj^'(-ic,i 0 ) 

The second term varies as (s)‘^2 anc | -j s therefore predominant. 


[E] r = 0(s~ ‘/2) s -j, 


[Ed has been computed to match the corresponding condition. We 
obtain from Equations (80) and (81): 

»— E n , - ia n _ SA (-ic,0) 

iti * ■ J sint/2 *■ j b * j; ~ x ° n rW'mc.io) 

♦ < l > (» ln -a) S ;sS- iC - 0i ) cost (97) 

Z n=l ln Rin ("ic.iO) 
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Thus, as expected. 


[F 4 = o(s ]/z ) cos<f> sint/2 


(98) 


The total electric field obtained from Equations (84) and (85) 
satisfies the edge condition. We now need to check the components of 
the total magnetic field, which are expressed in terms of the functions 

^ e 0n^ anC ^ ^ e ln^ ‘ ^ ecause the greater complexity of these functions, 
we will only determine their order of variation. Since the order of the 
singularities of each component is not apparent directly from the defini¬ 
tion of the functions, we will write out the diverging terms for each 
component as function of n and £, and then we can write the components 
of the total magnetic field H. Since we will only be looking at singu¬ 
larities, we will not include in the following expressions, the summa¬ 
tions corresponding to the incident waves since they are well defined 
at the edge, that is to say, of order 1 or smaller. 

For the n-component of H, we have the following singularities at 
the edge: 


For n odd. 



cose 


4i S5 n (-ic,0) 

, A 2 0 UV~r. “ , 2 2,5/2 
ckd Rg n (-ic,iO) (n +? ) 


[n2-f,2 + o( n 4 >n 2^2 > ^4)] 




. i4 
ckd 2 
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COScj) 
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(99) 


[5 2 -n 2 +0(n 4 ,n 2 t Z ,C 4 )j 


From Equations (89) and (90) we find that: 


and 




t 




cost 


1 A/ 4 2 2 .4, n/ -1/2, 

.7 o c/y 0(t) ,n f, ) 0(S ) 

(/■+>■: )- j/L 


(ioo) 
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[Hi 


2 ck 2 d 1/2 


+ 6 


S Qn^ ic -°) 

n=l ° n Rqp ^ (- i c, i 0) 
x ^s" 3 / 2 +4 (1-^On^) sin 2 t/2 sint 

, s 1 n( -1 c »0) 

i, b T" R (4)/ 

"in 


z n=l "■ Ri: ; (-ic,iO) 


^-3/2 t 4 s V2N s . n 2 t/2 sfnt 


cos4>+0 (1) 

(105) 


1 i 

Again by definition of e x and e z , the coefficient of the improper 
singularity in [H] f , normal to the disk, disappears. Thus, 


[H], = 0(s~ 1 / 2 ) cost, sint 


(106) 


For the ^-component of H. the singularities of the vector wave 
function involved are of order They are, however, cancelled 

by the coefficients introduced to satisfy the edge condition on [E]^. 

For n odd, 


[«4 4> i i ,<".'.4') = > 4 ^ 7 ' - - si ' n * + on) 


ckcj2 Rg^ (-ic, i0) n 2 +C 2 


(107) 


• - ~V s1n * + 0<1) 

ln * ckd 2 Rj^(-ic.iO) n+r 


where 


n = 1 17 
*+' 2 ?Js 


cost/2 


From the definition of the scattered magnetic field and from 
Equation (107), we obtain for the total field: 
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2H 


ck 2 d 3 / 2 


■ ", S 0n<- ic ’ 0 > 

B x S a 0n “(4) 1 , . 
n= R 0n ; 


:-l/2 


0) 


. .1 r' . S 1„<- 1C '°> -M, 

+ L b in T4T 1 - s 

n= R in ( _ic ’ i0 ) 


cost /2 


cost /2 


sin<f>+0( 1) 
(108) 


By definition of 8 * an< l 6 z> the coefficient of s'^/^ disappears 
and finally 


[HL = 0(1) sin<j> 


(109) 


Therefore, the total magnetic field H satisfies the edge condition. 
We have thus verified that the scattered field obtained by Flammer by 
matching the 4 >-component of the electric field at the rim satisfies the 
whole edge condition. 

We summarize the different results: 


[E] n = 0(s 1/2 ), [E] c = 0(s' 1/2 ), CE]^ = 0(s 1/2 ) 
[H] r| = 0(s- 1 /2 )> [H] c = 0(s' 1 / 2 ), [H]^ = 0(1) 


The formal solution for the bistatic, normal incidence, perpendi 
cular polarization case satisfies the boundary conditions on the disk 
including the edge and is therefore proved. In the next section we 
will compute the scattered far-field for the above case. 


b. Bistatic scattered far-field 
in normal incidence 


The bistatic scattered electric far-field for normal incidence 
and perpendicular polarization can be obtained from Equation ( 88 ) by 
replacing the vector wave functions by their large radial argument 
approximations, listed in Equation (34). From Appendix A, we will use 
the approximation of the coordinate system for large 4 . 
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When o 1, 
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c; - k r 

1_ » 1 
rd r 


Thus the scattered field becomes: 


Ej = - 1T 
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( 112 ) 
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We can easily check that E^H 3 is parallel to e r since we are in 
the far field region. 

The electric field given by Equation (112) has been programmed, 
and the results are compared with those obtained from the program de¬ 
rived by Hodge[9] with Meixner's solution in a later section. We will 
now study FI aimer's approach to the arbitrary incidence, arbitrary polar¬ 
ization bistatic case. 


B. Arbitrary Incidence 

In the previous section the disk scattering for a normally inci¬ 
dent plane wave with perpendicular polarization has been solved. The 
general solution for arbitrary incidence can be attempted following the 
same approach. 

The general formula for the incident field, Equation (3), can be 
separated into two components 


E^ = E^(cosa F 1 " + sina E 11 ) 

where 

F 1 " = (cos8 0 a x - sin0 o a z )e"^ ‘ r 
F* 1 = a y e - iRl * r 


(114) 


Let us define E 3 ! as the scattered field corresponding to the in¬ 
cident field F 1a . The general scattered field will be obtained from 
and E 3 ' 1 as follows: 


E 3 = E 0 (cos« E 3 " + sina E Si ) 


(115) 


In.order tp follow the method of the previous sections, we will 
expand F 1 " and F 1 - 1 as sums of.two different expansions. We will Uise 
Equations (44) and (46) for F 11 , and Equations (47) and (49) for E 111 . 
The H-fields can be expressed with the corresponding spheroidal wave 
functions. 
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(117) 


where 3 Y + ^ = 1 > 3^! + 6 Z = 1. 

Unlike the case of normal incidence these expansions involve all 
the varues of m and, therefore, all the different ^.-dependences. The 
reflected fields are calculated using the same term by term matching 
technique as in the normal incidence case. From Equations (30) to (33), 
we obtain the following expansions. 
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II II 

For arbitrary values of 3^ and b£, the edge conditions is not sat¬ 
isfied. As in the normal incidence case,_we will try to satisfy the 
edge condition by simply satisfying the [E]^ condition at the edge. We 
will study the perpendicular polarization case in more detail. The 
method can be exactly paralleled for the parallel polarization case. 

We will now expand the 6-component of the total electric field in the 
vicinity of the edge for the perpendicular polarization case. We need 

the expansion of the summations on and in terms of n and t 

In the same way that we derived the expansions of [E 1 +E y "]a in the 
vicinity of the edge for the normal incidence, perpendicular polariza- 

_ X (1 ) 

tion case for plane wave expansions in terms of M PA , Equation (80), 

— _/i\ c un 

and M e fT /, Equation (81), we obtain here in the arbitrary incidence 
case: ,n 
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( 120 ) 


where 6 mn is the standard Kronecker symbol. 
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x The condition on [E ] +t ]. is thdt it be zero on the rim. and 
. z must then cancel the coefficient of the singularity .;/(r,2+r2) in 
Equation (120). If we reorder the coefficient of the singularity and 
equate it to zero, we obtain the following equation: 


x r 

n=m 

L 


!ni-lnL~I 5 . >Q I-. (1-io h )(l+d lm ) 

c coso n (4) , . . , ^ 


+ y< f m+ln^ ’*0^ ~Wln^ c 
n=m+2 c COS6 o ^{^(-ic.iO) 


i ? Ym n( ° o) W~ 1C ’ 0) 

z n=m c Sin0 o R m ! (-ic,iO) 


: osm: = 0 


( 121 ) 


Note that the coefficients and tt have been taken inside of the 
summation of index m. 

Because of the .^.-orthogonality, the coefficient of each factor 
cosm<{. must be equal to zero. Because of the relation k^+.-^l, we obtain 
for and the various following relationships that must be satisfied 
simultaneously. 

For m=0 
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( 122 ) 
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for m=l 




y. Y 0n (o o ) S 6n ( ' 1c > 0) 
n=l C COS6 o R^'HcIO) 

2 y' Y 1n (e o ) S 1n ( - ic ’°) 'l 
n»l c s1 " e o R<J>(-ic,iO)y 


= 1 



(123) 


In the limit when e Q tends to zero, these coefficients lead to 
those of the normal incidence case. 

For m>l 
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(124) 


In order to satisfy the [El condition at the edge, B* and b^ 

<j> X z 

must satisfy Equations (122), (123) and (124) for every value of m. 
According to our definition in Equation (116), B^ and 3* are independ¬ 
ent of m since they represent the relative weights of the two differ¬ 
ent incident plane wave expansions. To see if a pair (3^,6*) can 
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satisfy the [Ej^ condition at the edge, the different expressions of 

i^, for example, in Equations (122) to (124) should be computed for 

various values of m. However, since the vector wave functions of index 
i=l are not independent, the fact that the coefficients should take 
different values for each m is not a proof that the solution is invalid. 
If this happens, another approach introduced by Flammer[7] can be used. 

Instead of one pair (b^.B 1 ). Flammer uses a whole set of pairs (,6* ) 
i 1 x z x m z m 

where B and ft satisfy the Equations (122) to (124) corresponding to 

x m z m i 

the value of m. Note that while ft, will then be the coefficient of 


the summation over n of M 


2(1) ; 


in the incident field, g v is not the 


weight of the summation over n of M 


M X 0 ) 


but of other functions whose 


^.-dependence is due to one trinogometric function only per component, 

— x ("M 

instead of product of trigonometric functions in the case of Mg' . 

This is due to the reordering necessary to obtain Equation (121). The 

new functions and M e ^]^ are defined in Flammer's book[10]. 

Their m-index, m+1 and m-1 respectively, represent the order of their 
-dependence. In order to write the incident field used by Flammer, we 
introduce the relation 


for m>l 


for m=0 
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(125) 


The incident field is then equal to the following formula: 
► 

E fJO , «> Y (0 ) / -i \ 
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(126) 


where b J ‘ and [>, L satisfy the corresponding condition. Equations (122) 
to (124)”. ” 
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A term by term matching technique leads to a reflected field which 
satisfies the boundary conditions on the surface of the disk and the 
[E], condition at the edge. But, when these coefficients are utilized 

in Equation (126), it is not clear whether these expansions still repre¬ 
sent plane waves. This question remains to be answered. 

By introducing different coefficients for each value of m, we have 
modified the initial plane wave expansions. Since the vector wave func¬ 
tions of index 1 are not independent, the incident field of Equation 
(126) may still be a plane wave. This needs to be checked for every 
value of e 0 . Computations of Equation (126) on the surface of the disk 
should be made to verify whether or not Flammer's solution is valid for 
arbitrary incidence. In his paper[7], Flammer does not mention the 
problem. It is in no way proven and a computational verification, 
though not a proof, should be made in future work. 


C. Numerical Test of the Normal 
Incidence Case 


A computer program based on Equation (112) computes the scattered 
electric far field in the normal incidence, perpendicular polarization, 
bistatic case. We will successively introduce the variables computed 
by the program, present the computed data, and discuss the results. 

1. Computer program 

The program, presented in Appendix D, computes the value of E, 
where E is defined by: 


pS _ ^o exp(-ikr) p 
^ " k kr b 


(127) 


The functions and subroutines presented in Reference[9] were used 
for this purpose. The output data are, as in Hodge's program for 
Meixner's solution[9]: 

the magnitude and phase of the e- and <j>- components of E, 


the normalized crpss sections of the disk corresponding to the 

2 *) 


e- and ^-components of ^ as follows: 


°a^ 0 s ,<!) s^ I 4irrC 

' ira r-*» 


im 4 

-Hx> 


lull 


= 7 ira.'v* s >i z 
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where a = 9,<j>. 
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(128) 





The normalization is obtained by dividing the scattering cross 
section by the surface area of the disk. The magnitude of the cross 
section in the scattering direction defined by (e ,<j> ) is equal to 

^(v«sK(v*s>]’ /2 - 

2. Data 


The results obtained from Equation (112) will be compared to 
computations made with Meixner's solution. Meixner's far-field solu¬ 
tion has been studied extensively by Hodge[9], who wrote a computer 
program which calculates the electric far-field for arbitrary incidence 
and polarization. The results of Hodge's program show very favorable 
agreement with measurements, with calculations based on the small disk 
approximation, and with Geometrical Theory of Diffraction (GTD) compu¬ 
tations. We will accept that program as a comparison source for the 
results obtained from Flammer's solution. 

Various plots of both solutions have been made and we will use 
the following notation for the curves, except where otherwise stated: 

-continuous line for Flammer's solution 

- dashed line for Meixner's solution. 

In the cases where $ =0°, H-plane, or *=90°, E-plane, the nor¬ 
malized crqss section and the phase will be tnat of the non-zero com¬ 
ponent of E, ox when 4> s -0° or a@ when 4> s =90°, respectively. For other 
values of <j> s , the cross sections corresponding to the two components of 
E will be plotted. 


3. Results 

From the plots drawn, we can see generally good agreement between 
both theories. We notice the following tendencies: 

The concordance between the two solutions is very good for small 
disks - ka<2 -. In Figures 5.2 to 5.4, the low frequency ends of the 
curves match very well. Figures 5.5, 5.8, 5.11 and 5.13 show a very 
good agreement for both magnitude and phase in the case ka=2. For ka=4. 
Figures 5.6 and 5.14 show a reasonably good_agreement for the E-plane 
cross section and for the phase of E in the H-plane, but the values of 
the H-plane cross sections differ greatly as 0] approaches 90°, Figure 
5.9. ka=2 seems a reasonable limit for the low frequency region where 
the two solutions match well. 

As can be seen from Figures 5.2, 5.3 and more from Figure 5.4, 
the oscillations of the cross section as a function of ka do not match 
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very well. Actually, the tendency is the same but the minima and maxima 
are shifted as ka increases when compared with the GTD solution in the 
backscatter case, Meixner's solution gives a much better result for 
large disks. 

For a disk of given size, ka, the agreement of the E plane cross 
section is rather good even for large values of ka as shown in Figure 
5.5, 5.6 and 5.7. The greater difference between the two solutions in 
this E plane appears to be the value of the minima. 

On the other hand, the ^-components do not match very well for 
values of o greater than 20° for large values of ka. Figure 5^.8, for 
ka=2, shows s good agreement but. Figure 5.10, for ka=10, shows H-plane 
values obtained from Flammer's solution much larger than those obtained 
from Meixner's as e approaches 90°. The difference almost reaches a 
factor 2.5 (i.e., ^ s 4 dB) for 0^=90°. 

We have seen that this difference appears also in the IT-plane 
cross section for ka=4 and in the 4>-part of Figure 5.12. The good 
agreement for ka=2 makes less likely an error in the programming of 
Flammers solution, but nevertheless, that possibility cannot be rejected. 
When (f> =0° and o s =90°, [E]^ is proportional to ei- An error due to a 
truncation or a precision problem may appear in the computation of 

g* for large values of ka. In our program we used the same criterion 

to truncate the infinite summations as Hodge's[9], and, therefore, one 
would not expect that this should lead to a difference. 

The phases, Figures 5.13 to 5.18, seem in general on good agree¬ 
ment for all values of e s , except for an evident problem in Figures 
5.13, 5.14 and 5.18. The phase of the 3-component of the field obtained 
from Meixner's solution is consistently 180° out of phase with that 
obtained by Flammer. Thus there is apparently a sign problem in one of 
the components of Hodge's since he does not obtain the same phase for 
the backscatter case, e s =0°, for ^O 0 and <j> s =90° as can be seen by 
comparing Figures 5.13 and 5.15, for ka=2, and Figures 5.14 and 5.17 
for ka=10. We verify that the scattered electric field must have the 
same phase for <f> s =0° and 4> s =90° when 8 S =0°: 


.=0° E 5 = E s e = E s e 


h =90° E 5 = E 2 e = E s e 
p s o e g y 


(129) 


Flammer's solution, as calculated, satisfies this condition. 

From the data, we see that we have an overall correspondence be¬ 
tween Flammer's and Meixner's solutions for the bistatic normal 
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incidence, parallel polarization case. However, except in the case of 
small electric circumference, ka, the curves are never quite equal. 

The weakest point of the solution is apparently the behavior of the $- 
component of the scattered far-field in the H plane in the neighborhood 
of the plane of the disk. 


D Conclusion 


In this part, we have rederived Flammer's solution for the normal 
incidence, perpendicular polarization, bistatic case. We have shown 
that, with an appropriate choice of the two expansions of the incident 
plane wave, the condition on the ^-component of the total electric field 
at the edge is equivalent to the whole edge condition. This remarkable 
property will not appear in Meixner's solution as we will see in the 
next section. It is a characteristic of Flammer's approach. No answer 
'oS been given to the question of the validity of the general solution 
for arbitrary incidence and polarization. The computational checks 
explained in Section B should allow a better understanding of Flammer's 
approach. The numerical test of Flammer's electric scattered far- 
field has led to mixed conclusions. Despite an overall agreement with 
Meixner's solution, Flammer's fields, show some important deviations. 
These might be due to computer problems. 

We will now consider Meixner's considerably different approach 
to the problem of the scattering by a disk. 
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Figure 5.2. Backscatter normalized cross section as a 
function of ka. 
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Figure 5.4. Normalized cross section as function of ka 
for * s =90°; e s =45°. 
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Normalized cross section as function of e s 
for ka=4; <j> s =90° (E-plane). 
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Figure 5.13. Phase of Ej* as a function of e s for ka= 
<f> s =90°. 
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CHAPTER VI 
MEIXNER'S SOLUTION 


The validity of Meixner's solution has been established. Theo¬ 
retical and numerical checks have been conducted in recent years, and 
they have shown very good agreement between computations and experi¬ 
mental data. In this Chapter, we will not repeat Meixner's solution, 
since the derivation is available in the 1iterati^e. The reference 
work used for Meixner's solution is Hodge's version[8]. Here we will 
derive the expansions of the fields in terms of spheroidal vector wave 
functions and show a proof of the solution using only vector wave func¬ 
tions without refering to vector potentials. 


A. The Vector Potentials 


Meixner used Hertz vector potentials to solve the scattering pro¬ 
blem. In this section we will give the expressions of the components 
of those potentials in terms of the scalar wave functions. Equation (7) 
This formulation is much more convenient than Meixner's for the compu¬ 
tation of the fields. Meixner's solution computes the potential of the 
scattered field for arbitrary incidence and polarization. We will use 
here the notation defined in Chapter V-B, Equation (114), for the in¬ 
cident field. 


E 1 = E q (cosa E 1 "+ 


r i 1 \ 
sina E ) 


(130) 


The definitions of E 1 , Equation (114), in this work and in 
Hodge's work[8] have opposite signs. Onr definition was chosen to 
avoid a supplementary minus sign in the expansion of E 111 in vector 
wave functions in Flammer's solution. This will lead us to introduce 
a minus sign in front of the scattered vector potentials corresponding 
to parallel polarization. As in Flammer’s solution, we can obtain the 

-SI 


scattered field in any arbitrary case if we know E $nd E^ -1 , scattered 
fields corresponding to the incident fields E 1 " and E 1 - 1 , respectively. 
We will define the vector potentials in order to obtain ESi. directly. 


Meixner's solution involves three different potentials, u 
the potential of ^he incident plane wave, tt is the "reflected" 


-i 


is 
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i s ] 

potential defined so that n +n has its components tangential to the 
disk equal to zero. Because of the shape of ni, the components of the 
electric field corresponding tangential to the disk to tP+P 5 are then 

equal to zero. tt s 2 i s a second scattered potential defined so that 
+usl +tt s2 satisfies the edge condition derived by Meixner: 


ft <*p> ■ 0 ■ s <'»> at <’ 3, i 

, 

ir x and 7Ty finite at the edge. 

where ir = ir x cos4> + Try sin<j> . j, 

also leads to an electric field whose components tangential 
to the disk are equal to zero. We list here the vector potentials for 

—C M j 

E i as defined above, using Hodge's notation 



(132) 

(133) 

(134) 

(135) 

(136) 

(137) 


i 
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/ 


_ ( CO 

$2 o v .-m 

"s = ~2 1 1 

iy k m=0 


[ X m + l +(1 - 5 m,1 )X m-l] cos « 

r H I' . 

1 !u i J . 1 + (l + * i)^ 1 tI sin 

—— i. u m+l m,l' m-lj 


m<j- 


cos-o 


038) 


The different functions used in those potentials are: 

“■ , n R mn ^“i c »i 0 ) ia \ 


mn 


S n,n ( - 1 c ' cos % )S mn ( - 1c>n > 


(139) 


1 (f»n«c,o) J 


.n R^(-ic,iQ) , n \ 


n=m 


S ” n< ' tC) R mn !( - ic - i0) 


mn 


S m^ c ' 0)S mn { - Jic ' T ' ) 


(140) 


U m = 


2 i 


.m-1 


W m-1 +W m+1 


^Vo 7 ^- 1 + Vl 


, = 0 , m > 0 

m — 


m 

Lui; = 0 , m < 0 

' m 


(141) 


xi ■ 2 i ”- 1 y . 1 ~y i , x" . 0 , m > 0 
” Vi\*i " 


V x i ■ 0 • " i 0 


(142) 


I 
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1 

I 


i 

1 


i 


, m > 0 


H , c . , _ v‘ £. S mr ( - (c ’ cos6 o )S mn ( - ic '°> 


L W m (c,0 Q ) =0, m < 0 


(143) 


»■ ,n 

(*'"' <C ’ 0 '> ) = nl t R<«>(-1c.10) 

mn 


<J> m (c,o o ) = 0 , m < 0 


, m > 0 


(144) 


As previously mentioned, the sign of Equation (133), is the 
opposite of Hodge's result in order to match the different definitions 
of EP". However, the signs of the coefficients of and have not 
been modified in Equations (137) and (138) in order to keep the parallel¬ 
ism between the parallel and perpendicular cases. The sign difference 
will appear inside U*' and x" instead. We, however, have the same u" 

i mm m 

and X^ as in Hodge's results. This is due to a sign error in the deri¬ 
vation of u" in Hodge's solution and to the fact that the x" are all 

m 3 II . m —c? 

zero. In order to simplify IT and X, their coefficients in tt differ 

3 m m 3 

by the factors (c o /k cos0 Q ) and (e Q /k ) respectively from Hodge's ex¬ 
pression. We now can express these potentials in terms of the scalar 
spheroidal wave functions. 

From Equations (7), (42) and (139) we can express V m cos m 4 > as a 

summation of functions: 

e mn 


V os m “ mkzj j‘ 

Urn n=m 

In the same way, from Equations 


r'J.’hc.io) 

(7) and (140), we have for 


’I'el^ (n »4>) (145) 
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m 



046) 



• n 


y 1 

L M 


N 

n=m mn 


S 

mn 


(-ic,0) 




♦elJ’ln.C.f) 


By using these expansions in the expressions of the scattered 
potentials. Equations (133) to (138), we obtain the following scattered 
vector potentials. 


1. Perpendicular polarization 


Si _ 


f. a> 

-T l i' m 

k m=0 


.n 


^ X rn+r (1 ' 6 ml^ X rri-l^ 


x ^— S (-ic ,0) mn 

n=m N mn mn 


R in )( - 1c ' i0 > (4) 

R ( *> Mc.iO) '' 0 ” n 


(147) 


Si 


r. « 


< 2 Jonl T ™ < <> ) «W|.1 c.10) ^ 


(4) 

mn 


*7! i-” 

k m=0 




». -n R^P(-ic.iO) m 

* i l— s (-ic,0) -99^- 

n=m N mn mn R (4) (-ic,i0) 6mn 

mn 


(148) 


2. P arallel polarization 


s« 


k coso q m=0 n=m 


03 “ y (0 ) R^ ^ (-ic , i 0) /., 
v v ron 0 mn ' , (4) 

l l 7T2Tf——~T '^err 


+ -^o — y 

k coso Q m= 0 
’i -n 


R^Mc.lO) ' mn 
mn ' ' 


i‘ m [u; +1 -(i + 6 lm )u;_i] , 


t 


l N S mn^" 1c, °^ 
n=m mn 


R^ ^ (-ic ,i0) 

mn ' ’ • 


mn 


R^Mc.lO) 

mn ' 


(4) 

; mn 


(149) 
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s« 


= - - l [u 1 ' , + (1+6, )U" ,] - 

.2 J n I L m+1 lm' m-1 J 


k cose m=0 
o 


<» .n 


’C-ic.iO) 


L t 5tt- 


n=m mn 


r ;„ i-'c.'o) 


fo (4) 

°mn 


(150) 


Those formulas will allow us to compute the fields simply by 
taking their curls. This formulation is very convenient for this pur¬ 
pose since to each scalar wave function in a given component, x or y, 
of the potential corresponds a vector wave function of vector e v or e , 
respectively, when the curl of the vector potential is taken. The ^ 
expressions for the electric and magnetic field will then be obtained 
in a very straight forward manner. 


B. Meixner's Fields 

_ || - 

In this section, we will derive the expressions of and E/ 1 
as summations of vector wave function. The fields are computed from 
the Hertz vector potential using the following general relations: 


H = i M Vxtt 
E = — VxVXTT 


(151). 


Using Maxwell's equations, we also have the relation: 


H = j^—v*E (152) 

o 

where Z Q is the free space wave impedance. 

Equation (152) will allow us to compute the magnetic fields from the 
trie fields. The scattered magnetic fields will not be listed 
*•< .! j-,<- of the simplicity of the transformation. In order to compute 
•'< *><• i>Js. we will use the following identity: 


k N 


U) , 

me 


where < 


x or y . 


(153) 


The final form of the expansion of the fields will be obtained by re¬ 
placing the vector wave functions N e * and N y by their expressions 
+ __ o mn 5 mn 

in terms of N p ,, and N e ~, . This way the electric field will be 
gm + ln Q m-in 

expanded in terms of functions whose :-dependence is indicated by their 
m index. These transformations are given by Flammer[10]: 


r 7 * 0 ) - 


N e * (i) 
qHI+ 1 n 


^m-1 n 


|m+l n 


No j 1 )] , m>0 

„m-ln 


Y(i) = + 2N +(i) Nn y(i) - ti 

Sfln ^''01 n » U 


(154) 


We will now list E 51 and p\ obtained b y applying Equations (151) 
and (153) to the corresponding vector potential, respectively. Equations 
(147) and (148) and Equations (149) and (150). 


1 . Perpendicular polarization 

e^ x = . I y ? Y (o ) V ( ~ 1c,lQ) gy(4) 

k m=0 n=m mn 0 R (4) (-ic,i0) 6mn 

. mn 

_ m J 

+ r F i' m (x" ,-(1-6 ,)x A ,) x 

k L n ' m+1 ml• m-1 

m=0 

i n R mn^("iC’i°) - y(4i 

X > S (-ic.o) -Tyr - N ( , xl4) 

N m „ mn v d r* / ^ -in\ e mn 

n=m mn R mn \-ic,iO) 




n-m mn 


R mn^(" ic,i °) 
ic,0) ^- 

R ■ (-ic ,i0) 


78 









We notice that the first summation is a "reflected" 
in Chapter V. 


field as defined 


Using Equation (154) in the above formula and rearranging in 
terms of Xj, we obtain: 




1 


CO OQ . 


’-illy. 

m=0 n=m 


k A Y mn v o ' ZVH 


(e ) R » ,( - iC - t0> 




+(4) 


m+ln 




♦f I • 

m=0 


1 ? l!_c , , r R ™Hc,iO) 


n=m mn 


[X 1 , N + ( 4 ) - x 1 N _(4) 1 

” +1 °m+!n °"-l 1 


with X^ = 0 for m < 0. 


(156) 


2. Parallel polarization 


pSIl _ 


kcose 


' I F v m „(e , ^i } C-*c,iO) 


o m=0 n=m mn 0 ~[7f) 


x(4) 


R mn ''ic»i0) e " ,n 


1 


kcoso 


I i 

o m=0 


.-m 




, 0 ), 


“ x(4) 
mn 


ffiVicJO) emn 


- Tf- s mn (-k,o) Mc ’ <0) n y(4 > 

n=m mn R^Hc.iO) °mn 

mn ' 


(157) 
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In the same way, using Equation (154) in the above formula and 
rearranging in terms of U^, we obtain: 


cSH = 


_.l _ V 'f („ , rn ,. , jT+(4) 

kcos «o m=0 n=m T ™ 0 Ri 4) (-ic,i0) m e - 


mn 


"m+1 n 


* %-ln^ 


t ,_L_ J 


kcoso 


o m=0 


R i^(-ic.iO) 

n=m t S "" ( ' 1C ’ 0) R< 4 >(--ic,iO) 

mn v 


[U» m + ( 4 ) 


(4)- 


11 M . n+s \ n" n 
m+1 e m+ln l m Vl e m-ln J 


w"' .h U =0 for m < 0. 
m 


3. Conclusion 


For an incident, field of arbitrary incidence and polarization, 
the total scattered E-field can be expressed in terms of E^ 1 and E 311 
as shown in Equation (115): 


= E Q (cosa eT 511 + sin™ £ Si ) 
= E o (cosa H 511 + sina H 51 ') 


(159) 


where, from Equation (152), H 51 = (i/kZ Q ) VxE^t. 

In this section we have found the expression for the scattered 
electromagnetic field in terms of the spheroidal vector wave functions. 
This allows us to deal with the fields directly instead of the vector 
potentials, and it therefore offers the same possibilities as Flammer's 
solution for studying the near-field problem in a much simpler manner. 

In the normal incidence case, we now have two different formulas, 
Flammer's and Meixner's formal solutions, giving the, bistatic scattered 
field. A term by term comparison of the two solutions should be carried 
out to see how it is possible to match the results. A first step how¬ 
ever would be to compare the scattered electric fields in the far-field 
region. This would be much easier since the remaining components of 
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E 5 are, in each solution, equal to single summations of scalar spher¬ 
oidal functions over n. The task of reordering those far-field com¬ 
ponents as summations of the Legendre associated functions has not yet 
been completed. The Legendre associated functions are the natural 
functions since they are orthogonal in the far-fiela region and since 
the spheroidal coordinate system approaches the spherical coordinate 
system in this limit. This orthogonality will lead to an infinite 
number of relations between the expansion coefficients dP n (-ic) and 
values of the radial functions at c=0. The validity of these relations 
could easily be checked on the computer since the different functions 
are readily available in Hodge's program[9] and in the one used in 
Chapter V of this work. A study of those relations as function of the 
variable c=(kd/2) might explain the discrepancy obtained in Figure 5.1 
for the backscattering cross-section. It would determine whether the 
differences are due to some inaccuracy in the computations of the coef¬ 
ficients or to improper truncations of the infinite summation^ in one 
of the solutions. The far-field of Meixner’s solution can easily be 
computed by inserting the large argument approximations of the spher¬ 
oidal wave functions. Equation (34), in the formulas for these fields. 
This showed a sign error in the 0-component of the scattered field in 
the perpendicular polarization case as computed by Hodge[9]. A sign 
error also appears in the e-component of the scattered field in the 
parallel polarization case in Hodge's work. 

In the next section, we will derive Meixner's solution with 
another method using vector wave functions exclusively. The edge 
condition will involve all the components of the electric field unlike 
Flammer's solution where the ^-component only was needed. 


C. Another Proof of Meixner's Solution 


As can be seen from the previous section, Meixner develops a 
solution to the scattering problem by using only one expansion for 
the incident plane wave in each polarization case - parallel or per¬ 
pendicular. His method is equivalent to adding to the incident and 
reflected field a summation of vector wave functions whose tangential 
components are zero on the surface of the disk but which has a sin¬ 
gularity at the edge. The variable coefficients of that summation 
are then adjusted to cancel the singularities in the ^-component of 
the total electric and magnetic fields, and to insure a proper be¬ 
havior of their other components near the edge. This is exactly the 
method that we are going to use to prove Meixner's solution. 

In this section we will first derive some useful vector wave 
functions expansions which are equal to zero on the surface of the 
disk. We will summarize the power series approximations of their 
components in terms of n and S in the vicinity of the edge. We will 
then be able to prove Meixner's solution, using the above expansions 
to match the edge condition. 
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I. V ector wave expansion 


This section is mainly based on the properties of the N x ex¬ 
pansion of a plane wave. Equation (47), for an angle of inci§?Rce 
o equal to 90°. 

a. Derivation of the vectorial relations 

—xM 1 

We will first establish some dependence relations for the N ' ' 

e mn 

vector wave function and use the ^-orthogonality to obtain relations 
for each value of m. We will then show that those relations apply 
equally to the odd vector wave functions. We can then find vector 
wave expansions whose tangential components are equal to zero on the 
surface of the disk and that behave like outward travelling waves 
in the far-field region by using a term by term matching technique 
identical to that used for the calculation of the reflected field 
in Flammer's solution. 

For o = 90°, Equation (47) multiplied by cose can be written 
as 0 

r v , (-/?) N x ^ 1) f'./,') = '6 1 1601 

p»-; ' a y 

m n r i r " nri 


From Equation M?) and Equation (14), we have, for n-m odd, 
< m ("/2)=d and, therefore, in the previous summation only the terms 
for n-m even will remain. 


) 

m=0 


1 ' W V2) N e (1)(n ’ C ’ 0) 


n=m 


'mn 



( 161 ) 


We will introduce here the functions N. 
a convenient ^-dependence, Equation (154). 


u (4) 


mn 


and N 


TT-(4) 


mn 


which have 


Equation (161) can be rewritten as follows 


oo ( 

l_ Y ln (v2) n; (1) 


n=l 


on 


f * i F ( ,+ V , .-i, n ( n/2 > "1 (1) 


/, i. 

m=l [n=m-l 


'mn 


♦ r 

i-wH 


Vi.nC/n « e (,, f -o' 

’ mn 
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(162) 


-< 


•r 


l. 


4 *» 


♦ 






where corresponding components of and N q ^ have the same 

dependence. 


mn 


mn 


A difficulty encountered with the use of this notation is the 
difference in the parity of n-m between Equation (161) and (162). 

Equation (161) is a summation of N*' ' for n-m even, while Equation 

_ f. ( i \ i \ 

(162) is a summation of N g k and N~ v for n-m odd. In order to 

avoid a possible confusion^the inde$ n p will be used instead of n 

in N*^and N'^ during the rest of this section. With that notation, 
8™ e mn 

the <i>~orthogonal ity applied to Equation (162) leads for every value 
of m to the relation: 


i (1+6, ){1-6 )y , (n/2)N 
=m-l lm onv m-lp 1 ' ' i 


p=m-l 


;(D 

"mp 


p=m+l 


Y m+lp 


(V2)N (1) « 0 
mp 


(163) 


Equation (163) means that the coefficient of cos in the n- 
and ^-components of the above summation and the coefficient of sin 
m$, except for m=0, in its ^-component are equal to zero. 

We can then substitute, for m greater than zero, sin mo for cos 
m^ and cos for -sin m<f> without changing the value of the summation. 
This transformation leads to the following equation for the odd func¬ 
tions: 





0 ) 

mp 


+ 


p-ii Tm+1 p 


(./2)N (1) ■ ff 

mp 


(164) 


for every value of m£l and (n»5»4>)* Note that Equations (163) and 

(164) are also valid for This is proved by simply taking 

gmp 

the curl of these equations. In order to match the behavior of out- 
wave travelling waves at infinity and still keep the behavior of the 
tangential components at f=0, we will make a term by term transformation. 


83 




I 


r r om Equation (164) and Equation (30) to f331, we ran see that, on 
the surface of the disk, the n- and ^--components of IT*'’] or fT^’l 


are proportional to R 
will therefore make 


R i- 1 (‘ 1C » io )» or R m+j^" ic ’ io ) respectively^ "We 
trie following substitution in Equation (164): 


. 0 ), 


3 m P. 


Smp, 


N +(1) 

§mp 


rw 

gmp 


R i-! P , - ic 'i° ) 

«Sp ( - (c - io > 


n +(4) 

gmp 


R (1 ! (-ic,io) ... 

"t+{j>: r (4) 

R m!lp ( - ic ’ i0 ^ 8mp 


(165) 


We will denote by F ( n, s, 0 ) the value of the transformed sum¬ 
mation. From Equations m 64) and (165), we have for every value of 
n: 


( 1 ) . . . . 

V*^ ■ pi, 

m-lp v * ' 

R (l) 


- X 


_+m “1 R ilL(-ic, io) 

« V 4 (n.s.*> + l /2) -rn 1 - 

&"P P-»1 ,n *" , R< 4 J-Ic.io) 


(41 


F satisfies for j n I - 0: 

gm 


C^ e ] (n,O,0) = 0 
0 

where ct = n ,0 . (167) 

[F ] (n,O,0) = 0 for rail, 
gm a 


In contrast to Equation (161) and (163) these F functions 

gm 

are not zero everywhere. Note that F is the reflected field cor- 

gm 

. .-sponding to the zero incident field of Equation (164). Since the 
•fleeted field is not zero, we immediately realize that the reflected 
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field is not equal to the scattered field in that case. We will need 
the expansion of the components of F in terms of o and e near the 

§m 

edge. In the following subsection, we list the expansion of the nec¬ 
essary vector wave functions. 

b. N^^andN"^ near the edge 
omn 5m-1 n 


We do not need very complete expansions of the vector wave func¬ 
tions in the neighborhood of the edge since we are only looking for the 
terms that do not have the behavior required by the edge condition 
for each component. The other terms will not lead to any problem 
in satisfying the edge condition. We will first consider the case 

°f , component by component. N^^has no singularity. Since 

S m+1n §mn 

Equation (161) involves the functions 'for n-m even only, we will 

_ + (-j) _ r ,• \ mn 

need N v ' and N ' 1 for n-m even only for the computation of 
8m+ln 8m-In 

F §m* 

Fo£ the n-component, we only need the terms of singularity greater 
than s~ a , where s is the distance from the edge as defined in Equation 
(89). 


n-m even,|n| << 1 and c£ « 1 


[n! (4 ! ] 


4i 


S mn^ 1c ’ 0 ^ nC {cos, 11 1 rtO( 

(168) 


gm+ln n ckd 


For the e-component of , we only need the terms of singu¬ 

larity greater than s~ 2 o m+ ^ n 


ity greater 
n-n even Jn| << 1 and c? << 1 


f, 


1 


rii+(i) i _ 2i S mn^" ic ’ 0 ^ e 2 -n 2 icos, , w <. . 1 *, 

Vm ^ckd 2 OT^io) } ] ° (S } 

" (169) 

For the ^-component, we need the terms up to the order of s* 5 
to test the edge condition on the ^-component. 
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F or '-m even, r ; 


1, c 


1 



Z. Sffln- ^ c »°) J sin(m+l) c ] 

kd 2 n 2 -K 2 |-cos(m+l) oj x 


x |( n 2 +-; 2 )[-(m+l)> mn +m(m+l)-c 2 ] R^J-ic.io) 

+ [l+m+ n 2 [c 2 + ~ (m+1) - £ (m+l)(m+2)] 

+ r , 2 [j (m-3)(A mn -m 2 ) -2m+3 J]— r j- T ’-1+ Q( s ) (170) 

c mn ( _ic » io )J 

Wp will need the expansion of ], for n-m odd 

5 >n+l,o* 


For n-m odd 


[N* (1) ] = nC 

gm+ln 4> kd 2 


4-A +m 

mn 


f ( A mn -ni ~ 2 ) +m ? 


sin(m+l )<*>! 
cos(m+l )4> j 


x s ^ n (-ic,o)R^ (-ic,io)+0(s 3/2 ) 


(171) 


It is well behaved at the edge since it goes to zero as s does. The 

case of N~^ can easily be. derived from that of 

gm-ln gm+ln gm-ln 

has no singularity. The n- and S-components of TT^ can be obtained 

gm-ln 

from Equations (168) and (169) by replacing m+1 by m-1 in the -de¬ 
pendence. The expansions of the '('-component of N" are obtained 

gm-ln 

from Equation (170) and (171) by replacing m by -m except in the in¬ 
dexes, sinfm+l) 1 . 1 by -sin(m-1)and -cos(m+l)4> by cos(m-l)<>. 


In order to compute the complete electromagnetic field, we will 
also need the functions M+M) 

and M~^ as defined below: 
gm+ln gm-ln 





As seen from their expressions in FI ampler's book[10], the n- 
and ^-cgmponents of these functions have a singularity of order at 
most s -2 and therefore satisfy the edge condition. The expansion 

of <(>-component of M 4 ^ is listed below for n-m even, 
gm+ln 


n-m even 


rfl + (4) t _ J_ 5 mn^ 1C, °^ -f c ? s (m+l)<S> 

[ gm+ln kd n 2 +C 2 [ sin 

!n{ (nVXm 2 -^) — R m )( - 1c > ,0) 


1 . 2 W m , r 2 X mn" m 

i + n —o- + s - 


2 n 


n \'] } ; f+0(s) 

cRi„(-ic,io) 

J 


mn' 


(173) 


We will also need t^ (1) 3* for n-m odd: 

gm+ln 41 


n-m odd 


L = hi S mn(" iC »°) R mn( _iC » i0 ) "I c?n( m+1 ) ' 

gm+ln 4> kd mn' mn' ] sin' ' 


f = 0(s) (174) 


The corresponding expansions for TT are obtained by replacing 

gm-ln 

m+1 by m-1 in the trigonometric functions. We note that the components 
of the functions of index 4 listed in Equations (168) to (170) and 
(173) have singularities that do not satisfy the edge condition, but 
that all the components of the vector wave functions of index 1 do 
satisfy it. 


F 4 is expanded in terms of N_^instead of N 4 ^ or N~^ . 

8m gmp gm+ln gm-ln 

From N 4 ^ , n-m even or odd, we can derive the expression of N 4 ^ 
Bm+ln gmp 

for p-m odd or even, respectively, by replacing m by in-1 even in the 
indexes and n by p. From TT' n-m even or odd, we can derive the 

-m 8m_1 

expansion of N ' 'for p-m odd or even, respectively, by replacing 
gmp 

m by m+1 and n by p. 
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F at the edge 

8 tn 


Here we will compute the coefficients of the improper singularities 
he components of F and 1/kv x F . The coefficient of 

r K ; Bm 

--i--F 7 o in [F ] is, from Equation (168) and (166): 

§m n 

/ v' M + , k, t , _ m ~ . ^ P 


- 

41 I „ • 

l (l4 

ckd lp=m-l 

v 


'ImX'^ora 1 


, , , ,,, Vlp M c ’ o) 

°m Yn, -'P R^f(-ic.lo) 

m-lp v ’ 


r' ( S m+lp^" ic ’ 0 ^ 1/cos 

E, Y m+1p (7l/2) 777),”". ."H sin 
p=m + l R.ti+ 1 p^ ' 1 c ’ 10 ) i ' 


cos \ 
sin m6 f 


. ?_ , 2 

fficient of — 15 ha ^ - 
(nS^) b/ ^ §m ’ 

wo components will, therefore, enforce the same condition. 

ifficient of in [F ] is, from Equation (166) and (1 

n + C o m ; 

c / 1 r 


of equation ( 


(175) 


(175). 


.» 

ent of -J * ~n in [F ] is, from Equation (166) and (170): 

rf+ C O m ? 

[ f p ( -n) u! p ( ;;- 

[p=m-1 R m-lp ( 1C ’ 10;J ' J 

“• , S m+lp^' ic, °^ , f-sin 

P Xi ’” + ’p (,/?) ^-TcTo, ! ' m) 1 «» ( (,76) 

be rewritten as follows: 

V S m-lp (_ic ’ o) 


(176) 


V , s m+lp^' ic ’ 0 ^ I'sin^l 

I.. w*«> **! 

tTH- ip > - 


(1771 
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This coefficient has a different ^dependence than the other two. 

The coefficient of _ ' ? Vx E 1, is obtained from Equation 

n +\ o m v 

(166), (172) and (173). It is equal to Equation (175) multiplied 
ky 2 * Equation (175) and (177) differ only by their ^-dependence 
and a factor m in Equation (177). 

Equation (175J_ and (177) give_the coefficients of the terms in 
the components of F and 1/k v x F that do not behave as required 
0^ o m 

by the edge condition in the vicinity of the edge. Note that for 
m=0, the coefficient in Equation (177) is always equal to zero while 
only the odd case of Equation (175) disappears. Thus F n and 1/k 
V . x F oo behave as required by the edge condition and, therefore, the 
missing relation in Equation (167) will not be needed since the func¬ 
tions TT and M are well behaved at the edge, 
op op 

We now have the necessary tools to solve the scattering problem. 


2. Solution 


The solution will be derived for the general case of arbitrary 
incidence and polarization. We will use the same notation as in part 
V for the electric and magnetic field. 

E 1 = E o ( COSa E 1 " + Sina E 1a ) (178) 

We expand E 1 and T’ from Equations (47) and (48) 


T ( " - 1_ I 1 r ro „<9o> N e 0) 

k cose 0 m=0 n=m mn 


r 1 * ■ A r y , (e ) jp'O) 

* m=0 nim mn 0 e - 


mn 


(179) 


From Equations (30) to (33), we see that, for n-m odd, the in¬ 
cident vector wave functions have null taggential components on the 
disk. Therefore the reflected fields are: 


r 


i 


R mn^~ ic ’ io ) -Wll 
Y (e J -r!n- N X(4) 


k c°se„ Lf, L mn'*o' „14), . . . ”e 

o m=0 n=m R' (-ic,io mn 

mn 


mn' 

to 


oo oo 


= l V 


mn ' ' ^(4) 


k m=0 n=m '* < "° ) rIJ>(-1c,1o) 


“mn 


(180) 
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Once again the "reflected" fields as defined here are not equ i. 
to the scattered fields since they do not satisfy the edge condition. 
We can now add to the reflected field any summation of functions F 
and still have a reflected field. The coefficients of that sum- 8 m 
mation will be adjusted so that the fields satisfy the edge condition. 


_ x f i 1 

From Equation (154), we see that N*' ' involves only even func¬ 
tions N + ^ and N~^ while 1 ^ invokes only odd functions 
e m+1n e m-1n e mn 


N* (l) and . 

°m+ln °m+ln 


Therefore to match the edge condition, we need 

satisfies 


only add a summation of F to F r and F to F r 1 sinrp r 

e m om c » since r 

gm 

Tne scattered 


the boundary condition on the surface of the disk. 


; e ■ l w ' 


o e: 


E S " 


-r" 


i 


k coso. 


V 

m=0 



(181) 


= E' 


k J, X F 

m=l mv o 


m 

where U mv and X mv are adjusted to satisfy the edge condition. 

The summation over m in F Sl begins with m=l since F is well 

°o 

behaved at the edge and does not correct the behavior of F r t We will 
successively consider the perpendicular polarization and parallel 
polarization cases. 


a. Perpendicular polarization 

Equations M 7 9) and (1801 can be rewritten as follows from Equ 
♦ io,i ! 541: ’ ^ 
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<T, f CO 

' ‘ Tc m^D | p=m-l 1 _<S om^m-lp^ e o^ N o n 


' £ +1 W 6 o )N o f 

p=m+l K mp J 

i v J v' R^l Me,io) 

? rA = .1 i] i (i+6. )(i-<s h .(ej-sJa.—:— 

k m=0 [p=m-l [ v lm M om' m-lp v o' R^(-ic.io) °mp 


I 

p=m+l 


■ R m+ n ( ‘ ic ’ io) —(41 

+1 Vip rm^T— t: "Cf 

+1 H R m+lp'" 1c » 10 ' mp 


(182) 


We must adjust the coefficient X so that the coefficients of 
the ill-behaved terms in each component disappear at the edge. 

Taking advantage of the <(>-orthogonality, we can write a specific 
condition for each value of m: 

for every m 10, atn =£= 0: 


k p=m-l ^ 1+6 lm^ 1-<S om^ Y m-lp^ e o^ N 

-r L ^p (0 °> *£’ 

p=m+l r mp 


mM) 


iilL(-ic,io) , 


+ 1 l y , (e ) #) 

* p=m+l m+lp 0 Rm+lp(“ ic,1 '°) mp 


(183) 


. R m- D (_ic ’ i0) —(4) 

' * p=Ll [ (W l« ,(M o. ) W e o ) rw’hc.Io, "“mp 


+ 1 X f (1- 6 ) must satisfy the edge condition. 
k mv om ' om' 
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We have seen that the functions of index i and behave as 

°op 

required by the edge condition. The absence of coefficient X is 
therefore not a problem. For any other value m, m>1, the con8¥tions 
on from the three components 'vV of Equation“M83) are identical 
.me" Me coefficients of the improper singularities in Equation (1681 
to (I 1 1 ) are proportional apart from 'heir ^-dependence. From Equation 
(173;, we see that this also applies for the o-component of the curl 
of Equation (183) which corresponds to the m-part of the expansion 
of the magnetic field. 


X mv satisfies, for mil: 


1 W-^q) 

<JO 

* pX, WV 1 

f 


p=m_l ^ + "i™^ 1 ““oin^m-lp^ <, o^ c 


R Sp ( - ic * io ) 


(184) 


+ X mv | pj’i (^iJ(l-6 om ) Vlp (,/2) 2 


p -■ n+ i 


c m+lp 


R m + lp MC ’ i0 ' 


R m .lp(- ic »i°) 

= 0 


We notice that, from Equation (42), we have, m ._0: 

i (e ) = 4 1 — S m , (-ic.cose )(l-s ) 
v om 7 lm 7 m-Ip v o N , m-ln * o om 
r m-lp K 

'm+lp< 6 o> ' fl£— Vlpl-'ncosep 


(185) 


Introducing the functions W and as in Meixner's solution, 
Equations (1431 and (144), we have for x : 


mv 


W m-1 " W m+1 
^m-l + ^m+1 


for m > 1 


(186) 


Witn those coefficients we can compute the scattered field of 
f i" pivpen 1 1 < j1ar polar i/at ion case. We obtain from Equations (180), 
f 18l) and ( H6): 
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E s ‘ - I ? ? (, ) R ^l ( ' iC ' ,0> *y(4) 

k m =0 n=m Y, " n 0 rFmTjo) e m„ 


Ri^C-icJo) 


. 4 ” L r }L ' , < rn) mn ^~ 1C> 10 ' tj(4) 

k L X m+lv N W~' c,o) R (4), ~TT N 0 
m=u n=m mn li' (-ic,io) m+ln 


R mn^(~i c »i 0 ) ' 

+ { 1 " 6 1 m )( 1 " 6 om )X m-lv l N~ W" ic,o) R 'f4')7 ~ ~T Vln 

n=m mn R' l-ic.io) 

mn 1 

(187) 

We can check this solution by comparing it with Equation (156). 

We compare the coefficients of and N~^ in both solutions. 


*_ . „ 2 j-m _ 2i“ m ,,,+m+l-l W rrf W m+2 _ 4 v 

r or m > u t ~ i X , — ■ 2 1 , , , 17 X ■» 

— k m+l k ^m + ^m+2 * m+lv 


*_ o 2 ,-m V A _ 2i' m 0 .-m-l-l W m-2‘ W m 

for m > 2 7 - 1 X . - —j— 2i -- - — 

- k m -1 k 0 +ii< 

m -2 m 


(188) 


k m -1 v 


The scattered field obtained by this method is identical to the 
one obtained by Meixner, as expected. 

b. Parallel polarization 

Equations (179) and (180) can be rewritten as follows, from E- 
quation (154): 

* - 1 - L{ i , (’n, m )(l- S Jy m ., p (e 0 ) h£> 
k cose m =0 ( i p=m-l r mp 


k cose Q m =0 ( i p=m-l 

+ P Xi Ym+i p {6o) "C' 


(189) 


r" 


k cose o m =0 \ p=m -1 


L) F (^^(l-^m-lpfe,) 

=0 I p=m -1 mp 


* f. 7 4) 

p=m+l K mp 
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We must adjust the coefficients U in Equation (1ST) so that 
toe coefficients of the improper singuT^rities in each component of 
the fields disappear at the edge. 

Taking advantage of the ^-orthogonality, we can write a specific 
condition for each value of m: 


tor m 0 at 




mu 


ITcos," p X, '* 

+ k cose n -m + i Y m+lp^ 9 o^ N e 
o p=m+i r mp 


( 1 ) 

mp 


1 


k cos0 


o p 


(190) 

w+(4) 

7Tv " n 

mp 


1 


k cose 


I rk •!„ 1 V. > IU 

l )Y -.(e) -TTV^-- 

.4-1 lrn om m - lp 0 ?T(-ic,io 

m- Ip' 

^ip Hc - ,o ) r < 4) 


o p=m+l 


1 


' , m+ip (9 o ) rr)) p , . . v 

Cipi-’'' 10 ) 


mp 


-- U F 


k cose n mv em 


st satisfy the edge condition. 


The functions of index 1 behave as required by the edge condi¬ 
tion and will, therefore, not contribute to the computation of U . 

For m greater than zero, the conditions on U from the three com¬ 
ponents of Equation (190) and from those of its curl are equivalent 
since the coefficients of thq improper singularities in Equations (168) 
to (170) are proportional apart from their <j>-dependence. For m equal 
to zero, the n- and C-components of Equation (190) and the ^-component 
of its curl lead to equivalent conditions on (I . The 4>-component 
of Equation (190), however, contains no improper singularity since 


the transformation of the coefficient of the singularity of[< (4) 1 , 

__( 4 ) m-l,n 9 

Equation (170), into that of En 4 3 leads to a m factor in the latter. 

mp 9 

It can be seen in [F ] , Equation (177), for example. For m equal 

e mp 1,1 

to zero, tiie ,|,-component of Equation (190) does not impose any con¬ 
dition on U in order to satisfy the edge condition. The satisfaction 
of the three v other conditions will then be sufficient at the edge 
to obtain the scattered field. For any value of m, the condition 

on !i i 


i’t.' 
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,Jl U +8 l» ,(| -'>o«>W“o ) c g ( . 1c ;io) 

v , e n 1 S m-f1p ( ' ic,o) 

' p=L w ■>’ c ^jy-ic,io) 


( 191 ) 


•Hi. 

'.i. SSs - "" 


Introducing the functions W and <|< as in Meixner's solution. 
Equation (143) and (144), and using the properties of Y mn ( e 0 ). Equa¬ 


tion (185), we have for U mv : 

W, 

U = — 
ov ^ 

W , + W , 

U„. = JllU- m -tl. for m > 1 


mv ih , + iii , 
V-l v m+l 


(192) 


With those coefficients we can compute the scattered field of 
the parallel polarization case. We obtain from Equations (180), (181) 
and (192): 


fS" 


OO 00. 


._ 5 s - »x(4) 

k COS6 o mk r.4 *" 0 e mn 


k cose 


00 ,-n Rmli-iCplo) 

r r 1 _ / v mn' y 9 

L ,Lsr 


o m=0 n=m mn 


* tU "*' v "Vln * °' S ™ ,U '"-lv "vj 


(193) 


We can check this solution by comparing it with Equation (158). 
We compare the coefficients of t (4) and r (4) in both solutions. 


m+ln 


m-ln 
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_ 2ii.. u- - 4ii' 1 *1 = _4_ 

k coso o k coso ~ k cose u ov 

o OlO 

->--m 0 .-m.m+l-l„ W +W 

for rn 0 -r -U" . = — c — 1 - IP ■ J H't- 

- k cos<) m+1 k coso & +<b 

o o m m+2 


k coso m+lv 
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for m • 2 U" , 

-- k coso m-1 
o 


2i 2i m-2 m 


k cose 


l 'm-2 v m 


k cose ' U m-lv 
o 


(194) 

We see that the coefficients match exactly. The scattered field 
obtained by this method is exactly that of Meixner's solution. 

We have given a proof of Meixner's solution by dealing with electro¬ 
magnetic fields only. The scattered field has been determined so 
that it satisfies the boundary conditions on the surface of the disk 
and the edge condition. This proof gives more insight in the behavior 
of thp fields in the vicinity of the edge than Meixner's. It, however, 
requires the calculation of the expansions of the spheroidal vector 
wave functions in the vicinity of the edge, and its edge condition 
consists of six conditions whose compatibility must be checked. It 
is algebraically more complicated than Meixner's, whose edge condition 
consists in one equation only since Meixner mostly deals with scalar 
functions. In the next section, we will calculate the scattered field 
for- n incident f ield equal to a single vector wave function. 


Scattered field of 


jle vector wave function 


From the solution derived in the previous section, we can easily 
calculate the scattered field for an incident electric field equal 
to a single wave function. The purpose of this section is to show 
that this scattered field contains all the vector wave functions of 
index 4 which have the same index m in their t >-dependences as the 
incident one. Let us, for example, consider an incident field equal 

to. Adjusting Equations (179), (180) and (181) to our case, 
K e m+ln 
we have: 




r r = 


i jC( 4 > t lki f 

k R in^" ic » io ) e m+ln k e m+l 
ron 


(195) 


Following the same derivation as in the parallel polarization 
(.tso, we obtain, from Equation (190), the following condition on U , 
• lu>' ' * the edge condition. 
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i S (-ic,°) 
i mn 


mn 


+ U 


m+1 c 




096) 


Thus, 


m+1 


S (— i c, o) 
mn' 


R*+ , (-1c,1°)[* m <* Kv2 J 


(197) 


From Equation (195) and the definition of F g , Equation (166), we 

see that the obtained scattered field involvlPs^all the functions 

^(4) 

and N ^ for the given value of m. If we choose an in- 
gm+ln gm+ln 

cident field equal to 4- N*' , instead, the scattered field contains 

k e mn 

all the functions and fT^ for m given, from Equations (154) 

e m± 1 n e m± 1 n 

and (195). Whatever the vector wave function of vector e ore, 
or derived from these, Equation (154), the corresponding icatte^ed 
field is not equal to the reflected field and contains the other vec¬ 
tor wave functions with the same ^-dependence. 

D. Conclusion 

In this part, we have shown that the fields of Meixner's solu¬ 
tion can be expressed in terms of spheroidal vector wave functions. 

We can therefore calculate the fields everywhere directly. In par¬ 
ticular, this will allow us to compute the near-fields and the cur¬ 
rent distribution on the surface of the disk. Meixner s solution 
is formally as convenient as Flammer's in this region. However, for 
computational purposes, Meixner's solution sh ° u1d be retained since 
it only requires the values of the angular and radial spheroidal func¬ 
tions and their derivatives for n-m even. On the surface of the disk, 
we can see from Equations (23) and (24) that we only need the values 
of the radial functions at S=0, since the derivatives can be expresse 
from them. We have also shown how the scattered field of a single 
incident vector wave function involves all the vector wave functions 
of the same ^-dependence. This is a consequence of the fact that 
the vector wave functions cannot be considered to be modes because 
of the dependence relationships. 
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CHAPTER VII 
CONCLUSION 


Considerable knowledge has been obtained about scattering by 
a metallic disk in the last three decades. However, most of this 
work was done in the far-field region. 

Meixner's solution can be directly applied in this case because 
of its vector potential formulation. The purpose of this work was 
to establish a valid solution everywhere in space. 

Flammer's solution was considered first since it leads directly 
to formulations of the scattered electric and magnetic fields. The 
validity of this solution, however, still needed to be established. 

We rederived it for the normal incidence, perpendicular polarization 
bistatic case and compared numerical results to data obtained from 
Meixner's solution. The general agreement of both solutions showed 
that Flammer's solution is acceptable in this case. We, however, 
encountered problems in the derivation of Flammer's solution in the 
arbitrary incidence case and its validity is still questionable. 

Thus, Flammer's solution did not provide us the expressions that we 
were looking for. 

Meixner's approach to the problem leads to a solution that has 
been, in turn, expanded in spheroidal vector wave functions. Another 
proof has been given directly using dependence properties of the sphe¬ 
roidal vector wave functions. This formulation of the fields is ade¬ 
quate even for near field computations. For example, the computation 
of the surface currents on the disk should be made with Meixner's 
solution. 

In parallel to the solution of the scattering by a disk, we have 
shown some dependence properties of the spheroidal functions. The 
derivatives of the spheroidal angular functions are not independent 
as the original functions are. Some dependence relations of the sphe¬ 
roidal vector wave functions of index 1 have been derived, and it 
is therefore impossible to consider a single vector wave function 
as a mode in these cases. 

Some questions remain after this work. The formal agreement 
of the scattered far-field of Flammer's and Meixner's solution in 
the normal incidence, perpendicular polarization bistatic case should 



be numerically checked as explained in Section VIB. The validity 
of Flammer's arbitrary incidence solution should be tested as explained 
in Section VB. If the validity of Flammer's solution can be established, 
it should be interesting to compute the surface current distribution 
on the disk from Flammer's and Meixner's solution. This could be 
a good test since the general formal solutions are too complicated 
to be matched term by term. 
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APPENDIX A 

THE OBLATE SPHEROIDAL COORDINATE SYSTEM 



Figure A-l. The oblate coordinate system. 


The oblate spheroidal cc 'dinate system has the components ( n , 
) whose characteristics are: 

- the ^-constant surfaces are ellipsoids of revolution (around 
the z axis, in our case). The intersection of those surfaces with 
a plane containing the z-axis are ellipses whose foci are the edge 
points of the circular disk in this plane. 
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- the '.-constant surfaces are hyperboloids of re.olution around 
the z axis. In the same way, the foci of the intersections with a 
radial plane are the edge points of the circular disk. 

- : is the usual azimuthal coordinate as in the cylindrical case. 

■',=0 represents the surface of the disk of radius d/2, c=h kd. 

varies from 0 to infinity whereas n varies from -1 to 1. 

The following formulas have been rederived from the general or¬ 
thogonal coordinates systems transformations as explained, for example, 
by St ratton[11]. 

We first establish the relations between the usual coordinates 
systems and the spheroidal coordinate system.. We then consider the 
limiting cases of large radial argument, c£, of the surface of the 
disk, £=0, and of the neighborhood of the edge, \=0 and n=0. 

1. Coordinate Transformations 


a. Components 
from Cartesian: 

f x = | \(]--r)(Ui 2 )] h cos $ 

> y = | [(1-r, 2 ) (l + c 2 )J % sin* (Al) 

i 2 = d r .. 

from cylindrical: 

. = | [n-n 2 )(l+; 2 )]^ (A2) 

from spherical: 

(A3) 

Sir... > I f 

1-n * / 

The metric-coefficients are computed from the following formula: 


= 0 2 +Z 2 = i 

... + 2 


cos.. = n 




r 






h. 

1 




+ 


2 2 

+ f—5■ ) where v^n, £, <P 

dv J \ 3 v J 


From Equation (Al), we obtain 


fn 2 + C 2 1 

. d 

rn2 + f2 1 

I 

OJ 

cr 

i 

^ " 2 

i- 

~+ 

spy 

PO 

i- 


[h* = | [(l-n 2 )(l+€ 2 )] % 


(A4) 


(A5) 


b. Unit basis vectors 


We give in Table A6 the dot products between the basis vectors 
of the spheroidal coordinate system and those of the common coordinate 
systems. The different formulas can be obtained using the general 
relationship in a basis (e..): 

t = l (a.e!.)e. (A7) 

The dot products are computed using the relation 


„ . ( 3a \ 1 

e a ‘ e 3 = \MJ F 


(A8) 


where a and 3 can oe x, y, z; p, r, 6, n, S, 4>- 
2. Large c Approximations 
a. Components 

Since Cis always positive Equations (Al) to (A3) reduce to: 

r 
X 

I* t 

cjl-n 2 (A9) 

r = f C 
n = cosO 

sine = J 1-n 2 for cS >;> 1 


= | fjl-n 2 cos4> 
= 4 f Jl-n 2 sin d 


2 
d 

p = 7 

d 
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This leads us back to the spherical coordinate system. The two 
main relationships are: 

n = cose and £ = (A10) 


b. Basis vectors 


For C£ >> 1 

U 


y 


-ne p + (1-n )e z - -e fl 


(All) 


\ - e r* 

On the Surface of the Disk 


The disk is represented by £ = 0. We have the relationships: 

p = r 


= 4Jl-n 2 


e + - 

n 


-A- A 

TnT & p 


(A12) 


4. Edge Coordinate System 


To determine how a component of the field behaves in the vicinity 
of the edge of the disk, the significant quantity of interest is the 
distance from the edge to the observation point. Let us call this 
quantity s in the coordinate system shown in Figure A-2. 


z 
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We will use the following notation: 


s distance from the edge to the observation point. 

t angle between the positive side (n>0) of the disk to the 
observation point. 

We restrict our domain to small distances from the edge: s << d. We 
notice that, in this system, with s <<d: 

t = 0 is the disk itself 


t = ir is the plane of the disk outside the disk. 

We will derive the relations between the spheroidal coordinate 
of small argument and this edge system. 

In the following, we have C << 1, jnj << 1. 

By definition, we have 


2 ( d Y" 2 

s = | p- 2 J + z 


(A13) 


In spheroidal coordinates, this equation can be written 
s ? = f j[(l-n 2 Hl-K 2 r s - if + nVj 
For n and f, small, we may write: 

fd-rfto + C 2 )]* = 1 - £- + f + 0(n 4 ,C 4 ,n 2 ,C 2 ) 


which leads to 

2 d 2 j ] , c 2 „ 2\2 2 A _ d 2 , 7 £,2 


(A14) 


(A15) 


s = 


i (£ " n ) +n = -yg (n +v ) 


Finally 


c d , 2 xf 2v 

s = ^ (n + f. ) 


(A16) 


<06 


* 



I 

I 

1 

I 

i 
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Gradient 

v <t> = 

Divergence 

v-F = 

Curl 

vxF = 


APPENDIX B 

DIFFERENTIAL OPERATORS IN THE OBLATE SPHEROIDAL 
COORDINATE SYSTEM 


1-n 2 

e + 1 

[in 2 . 1 

C 2 + n 2 

an ') d 

2 2 

C +n 

2 

1 


3 [(1 

-n 2 )(1 H 2 )] h 

at e t 




a- (L(i-n 2 )(«V>]‘* F n ) + Ij. 


Dn 

. 2.2 

C +n 


3_ 


2\ / 2,. 2\,JSr 


Jl-n 2 „ 




aF ' 

_ —1 

2 ^ 


(B2) 


IP 7 2 

1 


^ (J^U - ^ 


7(l-n 2 )(l+C 2 ) 


- at 


6n 


^=?-T7(F„)- ,===- - (Jl-n 2 F 6 ) 


7(l^ 2 )(1-n 2 ) 3 * 7T 2 ' 2 


/.v 


F'~ 2 M f/7*7 F c> - ^ ,+ 7 ^ 2+e? f ,> e f | 


* 

(B3) 
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APPENDIX C 

ERRATA IN FLAMMER'S WORKS 


The errata contained in this appendix have been collected by 
Dr. Hodge, Dr. Garbacz of the ElectroScience Laboratory and myself. 
The following corrections should be made. 

A. Errata in Flammer's Book [10] 


Page 42 - Equation (4-6-14) should be replaced by the following, for 
n-m odd: 


•n-m-1 JB . rn+1 ,mn , . , 

(1), i Tm\ c d, (-ic) 

V' -ic,io =- - - , 


(2m+3) l ' <"{-10 ifpU 


(Cl) 


r=l 


and Equation (4-6.15b) by, for n-m odd: » 

1 n - n,+, (2m + 3) l <"(-10 


R<^-ic,1o) ■ 1 


r=l 


cR^^-ic.io) 2 m m! c m+2 d n j n (-ic) 


(C2) 


Page 43 - ji should be replaced by 

n-m- l- j| - n R U) f or n _ m odd^ Equation (4.6.16b) 

Page 47 - The ^dependence term cos m( 4>- <t>*) should be introduced in 
the expansions of the Green's functions in Equation (5-2-11) and 
(5-2-12). 
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B. Errata in Flammer's Paper[7] 

Part I 

Page 1219 - Read i r+m_n instead of (-l) r+m_n in the definition of 
R mn^(' ic,i€ )» E< 1 uation O 4 )- 


Page 1221 - In the definition of a , 
side of the equation should be: 


2i m N 


on 


l 

r=0,1 



Equation (36), the right-hand 


(C3) 


Part II 


Page 1225 - The last inner bracket of Equation (11) should involve 
the derivation of the radial function and be written as follows: 


(D* 


In 


(-ic, i4) - 


R: (-ic,io) 

in R( 3 )(-ic,i£) 

In 


Rj^(-ic.io) 


(C4) 


( I \ * MV 

R on (“i c «i°) should be replaced by $^ n ' (-ic.o) as noted 


Page 1226 

by Flammer in the errata to his original report of this work. 


C. Errata in Flammer's Original Report [12] 


Page 12 - Replace sirs by cose in Equations (49). 




APPENDIX D 
COMPUTER PROGRAM 


This appendix presents the Fortran program used in Section VC 
of this work. This program has been written to compare Flammer's 
and Meixner's scattered bistatic far-field in the normal incidence 
perpendicular polarization case. It is based on Hodge's program[9] 
for Meixner's solution, which provides the structure and most of the 
subroutines used. The reader is referred to reference [ 90 for more 
complete details concerning the input and output, and for structural 
questions. First, the use of the program will be described. Some 
remarks will then be made on the programming, and a listing of the 
program will be provided. 

1. Use of the Program 

The program uses the same input procedures as Hodge's program 
except that the two following data, angle of incidence and polari¬ 
zation, have been suppressed since only the normal incidence perpen¬ 
dicular polarization case is considered. The inputs are made free- 
forrnat on a teletype as follows: 

1. KA = c, electrical circumference of the disk. 

2. THETA SCATTERED = e s [degrees]. 

3. PHI SCATTERED - ^ [degrees]. 

As in Hodge's program, any of these variables can be incremented. 

The corresponding inputs and the output formats are identical to Hodge's. 
The command "ESC" has been maintained. It interrupts the current 
calculation and requests a new set of data. No negative size of the 
disk should be entered and a disk of size 0 will terminate the program. 

Unlike the general case treated by Hodge, this expresison of 
the far-scattered field in the normal incidence case is here a double 
infinite summation over n and r, the index m taking only the values 
0 or 1. The summations are truncated in the same way as in Hodge’s 
work. The function involved in eacb-term of the summation are compared 
to a limit fixedgin the program, l(r u in this case, near the overflow 
level, i.e., l(r° for the computer of the ElectroScience Laboratory. 

To prevent any overflow, the summation is truncated when one of the 
tested functions passes this limit. This procedure allows a full 
use of the capacity of the computer and suppresses the need for a 
manual search for the optimum truncation. Flammer's solution involves 


(j 




the spheroidal functions for n-m even and odd and their derivatives. 

The subroutines computing the eigenvalue of the spheroidal scalar 
wave equation \ (-ic) and the expansion coefficients cr n (-ic) have 

been extended to n the n-m odd case from a previous work of Dr. Hodge. 

The subroutines computing the functions for n odd have been added. 

2. Remarks 

The following remarks refer to the program listed in Section 
3 of this Appendix. As in Hodge's work we will refer to a line of 
the program by LN. First, the different parts of the program are out¬ 
lined. The structure is similar to that of Hodge's program. 

LN 10-53 Input 

LN 54-75 Computations of the needed functions 

LN 76-101 Computations of the components of the E-field: 

- cross-section 

- phase 

LN 102-104 Output 

LN 105-108 Incrementation of the chosen variable. 

The needed subroutines are listed in LN 113-564. The notations 
of the spheroidal functions for n-m even as defined in Hodge's program 
[9] have been preserved. The derivatives of the functions are char¬ 
acterized by a letter " P' at the end of the name of the function while 
a suffix OD refers to the function for n-m odd. The following notations 
must be introduced: 

A and B are the two summations involved in the computation 
of and e*. 

i 

ALPHA and BETA in the main program are equal to B and 6 , re¬ 
spectively. Those names come from Flammer's definition { 7 ]. 

E0 and El are the summations in Equation (112) apart from the 
^-dependence and some trigonometric functions of 0 . The coef¬ 
ficients necessary for the computation of E0, and El, are com¬ 
puted in the and R(^(-ic,io) subroutines, respectively, 

under the name°C0F. n 


3. Listinc 


In this section is given the listing of the program described 
above. 








1 I'HTjOuS 3,t< 

.• r F Ai FlH.I S,rnTtHf'i BY a CinCUiAK FT T A(_ l IC ntRK 

3 f I (Ji # .-i/■ L T'lC I Ltf 'Jflt 

4 r I II H'.iilC'H h 1)1 AUlSATlON 

ilMfi mil r<s.si-.S»S 9 
t l l'-f'JUlu'j L7>' f L(b) ,VAR(3> 

7 CU'-PU '< F 4 ,F 4Ptf 0 .t 1 , I Y , a ,t *rOp ,F 1 HF ,f.PHI , Al.oHA , BE Tfl 

(• cU'’VI l t 1 (’ ( 3 'ill) »l < 1 00 .*>0 j rjFFb ( SO ) < F 4 ( 5 P > ,F4 p I FO ) .COF < bO > 

? 1. bn (f IJ) , S 0 | ( SO > «(M SO I ,pn I L>u ) . StT A ( tin ) iSI TAOO I 5(1) 

IL r>ATA L.'.HhL/<-.4F T HE. S PUT b / 

11 IAr(0..t.) 

12 VKTTt **%«■'(>) 

i 2 b F 0> • ■‘A I ()*«///«»< T Ypt I.f Sr» TO HFSTAPT PP.C6R Aw ) • t / , 

10 lMTYi’l. K/. = 0 TC ST UP PROC-OftM) • > 

I t C A| L k. C(V*»^) 

16 4 oor=i 

17 IU_l = ll 

In Ior>EX = t 

1° l Kt TF U*.«*7J 

,»0 ?7 F T< 1X«///,1*» * 1 . hA««l»y*'s ’> 

21 l<LAO(ti»-> vFP(l) 

?.?. JMtfArtl U.Lf.O.H-O TO 5 

23 UrtjTH <3 »«:<*> 

24 ?» FuuMi f(IX,*2. THETA SCATTERED = •> 

,'S FJ fc. ft & t o » - * VAF<2> 

2f- if.KTTM'NP'J) 

27 pb F Oti'iA r ( lX, * < , PFI SCATTfhED = *) 

2n REALM b*-) VARC3> 

2 C fkTTMo.oO) 

30 30 FOruiM U> »//«lX»*WHlCM V AKl APEl IS TO »F INCREMENTED?') 

31 RErO( b«-)NVAF< 

32 IF I «Mi/AR. L t .0 ) CiVAP.RT .6) >00 TO ?S 

33 WrtTTE<B,.H> 

-4 ? 1 FUR'aaT (IX, 'TYPE NUMBER Of CA<:Eb:»> 

3b PEaC(C'-) ROC 

WKtTK«A*J2) 

37 ?«i FO^'-irt F (IX , •UhaT IS THE I^C^EF'EIkt? ' ) 

3? REAi,MC.-M/1I'.'CRE 

f.V Wi\VAK = -'TV Ah 

4;, IL| l=2*Hl"VAH-l 

**i IUl2=TLI.l + l 

4/- wk r r e 1 3 tt-.u ior rflc ji n) »l<> bfcL (ill?) 

43 24 rai'i'iA r(i *,//.?x*2a-*lx«’ rKtsF ofctioiv* *?i x* »f mOrh»,/» 

44 113y ,« Sli', 1'1 r./ <pi*A» *> ) * , J.lv , *n E 1 A * « 15X « *i M t • « / * 13X « 

4 b 1'1 PF 1 A * • > X . *PnI * *OXt *:4AP* , bY . ’ PH A SE » . ‘ X . • ,-l AG • , 4X , 

4f- i 'Phase '♦ / J 

4 f b C=VAR(1) 

46 w hftAXs4S 

40 H|\J*. AX — 4'j 

SO jKI»v|AX = 4b 

M Tm- S = vH2I*a.lH159/lflO 

52 fI xsCgSF lntE» 

53 FHTSti/'M ( 4) »;x. 1.41bQ/lgo 

54 (>U 10 F4=l,2 

bb F'SHM-I 
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>j6 ru-o 

5? CA( l O n t ib|» ec ,w»Wf''»‘AX) 

St CAi.i u^rtf-i"(c A>.,r>N«.ax, r p f<MAX i 

V ' C - tHff.) 

faO bb CA( (_ 6 “Tjt P ( f , f f ’MA X ) 

ht call ol^aoi ( c, w , a , i'),mf ^a*.ph^ax,ifr^aX, 

b;-. Call k-. L*r (f t/\*«* Ink,-i A x) 

b/ CA ( L OPANGtL'tJpAX * 1PRIWAX) 

to« CAlL c 1 M-lili ( NP^AX * c o ) 

bt: GO TO "7 

bf bC. CAI L bl''MlMr.,|^.KAX ) 

b7 CAI l U L ’ E «- fJ ( C t M,N('"Vx) 

r,n Cam v 'hi r )<C ,r,3, H ,P'«VAx,i.M""AX t tRRma x ) 

bV Ca' L PMLT.'»CETa« •"» fbrti-'AX) 

7l) CAI L OiVN(, (f ftMAXtXPNF'iAX) 

71 CAI I PltLljJ (N^«AX,El) 

7i- t>7 CONTINUE 

75 1 r, C Of. TI riUfc 

7“ Ac»HAsiV(t3-A ) 

7 j Pt.TA = -' 1 /< Fl-A) 

7b f r 11 L = rt L P 11 A » I 0*SlMP.US> 

77 PhhI=ALphA4LG*PTA*C0S(PHts>-' 3E 'TA*EI*S i ^»THES>*''0s<PHT<;) 

7A p El-i*' GTsC Ai<S (E ThF' ) 

7? Em/>oP»CAPS< fF-Hl> 

Oil SIi;tiiL= ( Ei v / <*T*?/C) **<> 

£1 Sir,PHl = (t.i v A( P*?/C>«*c 

ap EK=REaL(UK> 

f l=AlMAG(tTHf) 

CM IFCfO.KP.G.) GO TO jfc 

fib AhArLI/Eb 

fib f'PP AT = 1 Aii/ 4 . xul 59* A TAN ( ApG ) 

»7 iFULK.Ll .( , ) .A'A'. l(. I.GT.u, I ) EPKrtT=FPGAT + 180 

»(■ IM (t K . L T . . ) .A ■«!•. • <tl*l T.o. ) )EHHATst PbAr-lflO 

fi9 lb F>»/lT = 9f) 

')•' IF f K'l •!-T.n . ) t PHATs-VC • 

91 l 7 EKrBF./iLlT.PpO 

9P F I r: A11’ 1 A G (t P F 1 1 ) 

9p IFrCP.PO.tJ.) GO TO IS 

9'» AKGsFI/FK 

?b EPpaF’ 3) f>0/3.11* 159*A TAfj ( Apt>) 

9b IF f f t'K.Ll ,C . > , ANp. <fTI .GT. u, ) > EpHAP=rPPAP*.180 

97 jFf (PF.lT .0. ) • ANI/* <£j .1 T.U. ) jEpnAP =r P , tAp»18l1 

9b GO TO 19 

99 j8 FPhaP=9i) 

1 .>«» iFMc-I .LT.O. )tPHAP = -90 

101 19 TKiJVrtP.tL.l'JA'MVAKsl 

lo* lkite< a**5) vr.ntW'VrK) iSicrHF.si&pul ,f i v 'A&r.ep«ar* 

It'-5 IFM/igPWPHAP 

10A FOr. vmT < jX.f 7.? , i* c ly, E It ,.* ) .2 ( 1 a ,E1U . 5«1 v ,F7. J»)) 

10b 7,7 IF ( INUt X*Eb.N0r>l’O TO 4 

10F J.\iPLx=irii’Ey*l 

10 7 VAP I IMI*V'AA ) S VA(? <NNVflb )*.VImcRF 

10b GO TO « 

109 5 CA( l. EX 11 

lip 4* CUwTlijUn 
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111 
11.' 

I I •> r 
3 i r 
U‘> r 
1 U r 
117 f 
11''. 
n< 

i^.t 

121 

i a > 
io> 
124 
i.*U 2 
13?., 1 
127 
1 2.0 

I. i r 

III 

I I, r 

111 r 
l 04 r 

IS-J r 

li'.! r 
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lit 

ll-j 

lHu 

14) 

142 

1**! 

lui» 

141 
14b 
14 7 
14c! 

1 4 9 
Il¬ 
l'll 
I‘j<? r 
IjJ c 

1 J4 f 

lit 

i-:.~ 

IV. 7 
lie 4 
l.i 9 
lt-ii 
lil 
1 ».» •* 

1..1 

lb4 

lr-i 


00 TO 4 
F Nil 

fill Alt. .s(h[F('TPAL MuuLf" C-UMCTif.A', S, OP AfcOmENT n» 
OK.ipfi u- t fVFMl or To CKUF4- = 

It'. jAL TO p"f (01) 

SUi ^out li.i (ii > 

CO-Plt.-V F'>,(Lp,C0F 

ro '’ y - n '* Flb( 1 r. o J »l,( i !)n ,50 j ,0r FG(so) «F4(5( ; ) ,F4p(«cO) ,CnF( *0) 

1 • ol> (:>U ) 

SO(1)=1 

TF JO 1 

nu 

SO ( 1 ) = ( 2*i"|- -1 )*$»< 1 ) 

rtl .< ».I. S J ,r.::;iAX 
l' = :>+ (iji .-A ) +■> 

V U ( l"'j+ l » — — | i. + m + I ) * $0 (l|F!) /i IV-M+k 1 

1 COI\| x!'JI IfT 

ft-THrt'l 

Finn 

OtiMVATIVf. (F ntll/.Tt SPHrhOlpAw aMOULaR 
FO-ir.TlOt i OP AHOO..EM 0; Ot'-^F.il f.',N, 

*** 1 Ti, rt-N OuOS 

SO., o JUT II, t Sr,A'OF , <,«i,i'jf. ( MAX) 

0 0- ,1,0'J FI ,(U#)|I'(1U0,'.0| ,C'NrG(sO) ,F 4 ( 5 r ) ,F 4 (»|K 0 ) ,CQF(S0) 
i.sr (4U) ,GOMbo> 

CO -PLI.Z Ml,FHpiCOF 
Sui (11=1 

MJr.vna 

(0 2 rtM = l,Nfl 

Soi (Us|c*|v>,.1)*S0t (1) 

9 CO’-rlUOF 

To j MN=1,NNi»aX 
f' , = 0*fJN + i J '"J 

SOp (H; j*l )=- ( A +fl + i.) »$c p (flf]) / (N"f.,*l) 

CO- T I i\lUFT 

He.TllRI'l 

CNp 

Obi aTE SbHFbOlRAL EIGcNVaUOES OF AKGlif'E.NT C, ORDER M,n 

SO, nOUTir e Ot EIt»MlC»Fi,NNMAX) 

C Or. '|0,’j Kif.-(ltn) 

0l"FMi,l0N XPibO 1 ,F ( 60) ,6t PbA f bo ) ,bLTA (50) 

co». ri'j"E 
r «c=.r >c 

7cr si.tiE-fb 
[Jill'. sfJNI'.AX+.J; 
i\ 1 

PI1)=1 
IH(1)=1 
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16 6 


iS-n 

167 


1 COhTliJ'JK 

lob 


no 2 1 ) 0 = 1 , "'N 2 

16',. 



170 


im=M 2 + 2 *iufe+is 

17 J 


Ix = 6'2 + 4 >* 1 , . + 2 * t s *1 

17c 


Alp* A < ! ia ) = (2*IV-i )+?*1 v* (IV-1 )-1 ) )/ t IX* (I x -*») ) 

17 A 


i- (•■ + x V- 'l t * ( ) V*'? • 

17h 

? 

r, t- T a ( i'jo+i ) pr?/Ti*^Q H T i r v *(xv+i(»ih*(iu-i)/(tx*ix- , *. 0 )i 

17‘w 


Pt-TA (iJi'!.-: + l )= 0 . 

i ;6 


no-A>r$( ALPf.oj i > >+r.f-s«tnT fl < 2 ) i 

177 


00 6 To'J=2,Vl!2 

175 


Atir ABS (f-L T A ( IOO) ) +1 B S t AI.ph A / r fit) )+APS(BrTA(XOfl + l> ) 

17. 


PLT/\ < 1 >' 9 ) rl f f /\ 1 XOO) *t‘CTA 1 100) 

If; 0 


IF | Ari.t'T.bO) B0 = AU 

letl 

X 

ooNTuiur. 

102 


AO = -ao 

105 


bot=bo 

lew 

15 

cumtimuf: 

Ifaj 


f’Ornoi 

186 


f )0 i;o lOosi , 2 JMI 8 AX 

187 


1 = 2 *IjO- 1 +IS 

186 


N =1 +M -1 

186 


a=ao 

ism 


B=po 

191 


ItWRr.l 

192 

?1 

I ls =0 

19? 


CU = (A*B)/2 

194 


XPfCU)t>o*22t50 

195 

50 

f Ko = (i)-A)/ABS»CO» 

19c 


ItPRsiKKK 

197 


IF ( It ) '»C ,'*1 * tt l 

196 

01 

; v lIr T f ( jf 

199 

02 

FOi.yiflTIxX, * 1 TfPAT1(i!vIS EXrtECFD PCR ElGtCVALUr ♦,X3> 

2 0 0 


60 TO 700 

201 

00 

IF ( r9K-ACC)24,29,22 

2 IV 

?2 

F>(P )=AL PG A ( 1 )-C0 

20 5 


TO s i= 8 *M 

2C8 


fl I )r< A, F-pA< 1-1 )-Co-BFTA < I-l)*lP(I-?l/F‘(T-l>n*P(I"tl 

205 


I'rtAOrA lSIPf I ) I 

2 06 


IF ( PMAi .,0 r. ) .0FF35 ) ^0 TO 7 

20 7 

5 

CO'lTINlJK 

20b 

J 2 

ooNTiiNf'F: 

209 


CO 6 1=2*NI 

21 C 


lF(fMl>)) | ‘,C,9 

211 

6 

TF 

212 

1** 

IP* I )=-l 

2 li 


GO TO 10 

21 *. 

o 

ini n=i 

215 

10 

IM1F’<I»-1F’(J-Il)*>.ll«fc 

216 

H 

I15= = I1S*1 

217 

6 

COntIM'E 

2 l*> 


IFIX1S.-X<.t>lfc,lb»l5 

219 


A s rn 

22 v 


GO TO 21 
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221 

1 fa 

f'=CO 

22 ; 


‘0 TO 21 

22 3 

?*» 

t’U = CO 

224 

700 

- 1 : . ( L ) = - 10 

*25 

20 

f Or. T IMJE 

22c 


JF(TS) rtUO.tOOiOUX 

227 

• 0 f| 

15 = 1 

2?e 


GO TO 1 

22V 

n i l 

CO'IT IIJ'JK 

2 jf 


K c. T UR ‘4 

23 j 

7 

M\lvAX=l-4 

23F 


f 1 1 = NNi"i AX* 5 

23* 


f-0 TO 4 

2 34 


3 Mi) 

23S, 

r 


C.V, 

r 

03| a Tt S^Hf'OinAL t I6E1 FiiWCTTOim fX^AMSIOii cOrFF IC I ENTS 

237 

r 

CF AKoOrf-NT C: OHUtH m.n.r. 

2 3:'. 

c 


23V 


$U»I.OUT U-t ChfPFMc.f iM4|F'AX,rM(iAX* iRpHAX 1 

24 0 


CO-.,»O.I tI(.(loo) «t'UUb,f.O) 

24 i 


r Wjvis ion iji (ioo) 

24 2 


C<i = c*C 

2 43 


i*i= i+i 

244 


A X^ = ? ♦! "A'F'Ax 

24b 


GO i wNsl ,Nf'|WAX2 

24& 


r. =c/,+ r-j rj - 1 

247 


J5= (-1 ) **IW>J 

24 b 


IK(TS.I-E.O) is=o 

249 


S = 1S 

250 

4 

l‘F>( T»K 'AX*3»sO 

251 


f>( imc- AX+2)=1 ,ut-?o 

2 52 


n ( 'u,, i j =0 

255 


l’ l f M * 2 ), = 1 

254 


v IO=<U'vi-l51/ <: + l 

255 


eo i.o7 Lu=i*ippnflx 

2 5c 


l .^1 L ”1 

257 


IF (1 L.or.oJ) l =IRR''1 AX+jj.li. 

2 5rj 


IH = 2*L'* - I6 

25 7 


lnrsy+lH 

266 


AR= ( W+J (J !-. + i.,*,l"i + lKv+i )*C?/| ( ?*iH l k+3,. t 2* I HM + 6 ) , 

26 l 


f'K- ( 2*X(>r , .» ( ) Rn + l) -2*, ) + c ?/( |2*irw.«l)* 

262 

j 

. t>* i^w+i))-Tr"’* i im "*■!) 

263 


CK-IK * 1 1 6-1 1 •+ r /■ < (2*IKM-.X)*|?»1RM» 1 ) j 

264 


1HLI* )j>H5»10b*106 

265 

1 0O 

0 (c. r i * L + 3 ) :-fCi;* , '(4' t o +J ) + ( Hr+Eir (Wf,) ) *t) f f<i.J »l +21 )/AR 

266 


niA:,=Ai:Pti’("'W,t *i) i 

267 


IF (i)f-'t7.6 ,41 . i. .OF + 3U | vjo TO 3 

2oJj 


60 TO 111 7 

2(,m 

106 

rF > fl+11=“<Ai-»PP(L*i) + (r , 4+t I6(7i|M))»0P(L + ?))/Cp 

27o 


cii/K,=iuis(i.i (u + 1 n 

271 


IF (OM.U-,.,1 .1 ,0f ♦5l'|fc,0 TO 3 

272 

1 0 7 

COnT I r;i i); 

275 


rL:MMil()t.'i v i 1 |tll » 

274 


t’U=ALoi-U>(l’L) 

275 


PLP=A’I j ( OF J OJ+11 ) 
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27t DUP = M-°&U‘(l.H.P» 

277 ''LsJIM'II-I 

t* 7 nct' r Ai.i 9 (1-M.P I 
279 rtrpl .♦‘"LP 

2Ao IF f i)l_ . <- r • 30 • )f>0 To b 

•L, i ro ; = i< f 'J«J+i> 

2dc = 

*Cw If- ( .'.C'J l! .l>r .1 .of +32)60 TO 'i 

2t-9 CO 11(3 = IPRMfcX 

2tV- 1 1 r j*2)=cr,M*P°< J+2) 

2P. 5 F = 1 

207 IS.c=Ii> + i'1 

200 in I'm) iyp*i9e»i y 9 

2A9 i 99 (.sb.o + M- lS* 

2V,; DU 1.1.0 1 = 1, iSM 

291 H' P=P*fISM-l) 


1 9». SU ' = 0 

293 P'|-|V = INRMfiX + l 

291. CO no 1 = 1, l '«IX 

29b lK=P*X+rs 

29f Su«=SU l! + f-*0<^N,I + l) 

29/ TK(I-J«I) nr- ,197,1)3 

290 197 PlUf’sP 

299 11A F={-F*(1K+2*M+1S-1)>/(JR»IS) 

300 ALCsPN f1 /Sl.‘fc 

3o 1 DO 119 1 = 1 ,1' 1 MX 

302 n=Ai_F*'.HNNiI + l) 

3;i A t Xu COi-TI'm'IE 
3,19 1 C Oi'TI im-je; 

3nb PtTuR.J 

3(ir> 3 1 KO,v,AX=LL-1 

307 <-0 TO 9 

301 ) ? IKBmAXsiKKpAX-I 

309 C-y TO 4 

31fl b NN-»'rtX=<M^-l )/J> 

311 HtTUKW 

312 turn 

313 r 

319 r f'tRAtlVC 0 CUFFFICTEImT StJhftOUTXwF 
31b r tVF.N. 

3lij C 

317 SUi(HOUTI lv t L’f'JtC C ,f**NMflflX ) 

3 If; tO^MO' 1 . f. l*( l0fl> » ul 2 OOifcO 1 ,UNrG(50> 

319 CO u i'm=l,Nr.,v,flX 

32c C2=C*C 

321 Mis2*-9N-1 

322 It- («,«f ,1> 00 TO ? 

323 OU b nM = 1,Mmiv,aX 

329 M'=2*l\l';-1 

3?b $ tNr(.(f9M»=U(''K,U 

32b f-0 TO 3 

327 a f’Arl.O 

32m 92=0.0 

329 Fi=riw«fJK| 

330 I'O 1 IRR=X,* 
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331 -? 

34,' Ah-ij* +i f +.'?)* f «s*“i+l 1 >4-l) *t 2/f <<£*,+2 *Tp(+S) 

4'; a* (;.*!+•'* t>\+-> > 

334 pK- ( r + | ’! ) ■* ( - ♦rR*3 1 -(•. 1 - ( ?* I r -*"fp ) » {,v4 J p + 1 ) -2*r-' ) 

i,i‘ i*cr>/< (,♦*:* i r-i > * ('•£**•*'*■* rh ♦3 > > 

4 36 CK- ( IK ) * ( 1‘i-i ) ‘C2 / ( ( iJ*[*2* j r. 4 ) * ( 2*y + 2*TR-l)) 

347 ^3 = 1'? 

330 t'2 = (U 

33y 1 C|= (BR + ri2-0M*n3)/Ap. 

34p A=n ( Uk> . 1)/81 

341 IWfc < f'jf 11 = A 

34; T uSitj' (L')U'(NM> > 

344 IF (rjilrt.l.T . 1 .1 F-3b)r,0 TO fe 

344 4 f 0|!TI.|O£ 

34t- A Rt-rijRl'J 

34m (■ l'!Uf 'lXr.NW-1 

34 7 At. TURN 

344 t'No 

34 r J r 

3SC r Otli ATt •jf'i'FXnOAC '<APIAL FUOrTIOf'. *14) OF Af.<;uwENT Cl 

4bl r CKOFH l- I TO fi-" fcVl.v- l.p TO DKLiFR l\,=>: + t*nWK H*.Z , 

3f)2 r A|_RO l,(,hM/> L 14AT14'i'i FUlvC T T <j ! V • M. 

3'ij C fL^O COFMirifNT CIF POP C OF d| JTaTION OF Fl 

3f>4 r 

od'j sU'.'»n»u’ I.M6. npf, au( c lr« “iF-Mfix,f.f.i”iA> • ipk-iaV» 

3bo CO* y,u>, f lu (10(1 > ♦!•* 1U0 , ) ,Of.F6(*SO) «F4 ( 5n I ,F4p (^0 ) ,C()F ( bO) 

337 1 tJ>>, (hO) .SCPtbO) 

3bl> CO-lPLt. v iX.F4,F4t£!»C0F«F4P 

35v lx~|P.0,1,0) 

3 tul FF/\C = 1 

36l r AC = l 

3m FAC^sl 

3bi C-K.. = l 

344 IHw.LO.O) f-0 TO 20 

3f>‘> r rty,V| = i;*l 

3fch 00 14 1« = 2,MAXM 

347 IKsMM-l 

364 i,rto = {F. *H -1 )*<2*J "I*OKO 

3 „•> FF<\c=i2>l*--I)»PFAC 

37c F Ac 2“ I ' »I M - l ) * ( 2♦ I *; ) *F AC? 

371 FA(; a IF>*fAc 

372 IF ,FF*Ar ,l,pr*l 7)GO TO 4 

373 IF ( r.",C.1.0F + 3U ) GO TO 4 

374 2r, |*=(0.0,0.0) 

37b 10 17 '4\l=l ,Nt;MAX 

37fc> FA=>*M"-1 

377 Ai?*|ij"'-i)+r 

376 S-U' = 0 

374 OKrliKO 

36'* F 140)61 = 0. 

301 Oo 14 ' f| = l.lRRM*X 

3n< IKr2*<'R”l) 

3a 3 r.y?p = (,c*i , »i,F ,nP > 

i«4 SO'fsSul' + SlIhP 

30 b I fiAbsnl’S t'Jl f'K # NKl I 



380 

a 7 

38 '1 p 
3(t l » 

39 t 

391 

392 i 6 

393 

* ni» 

3‘)s 

39* 

39/ 

396 
399 
4 01. 

HOI 
4 ,2 
4 J3 
*♦04 
Ht6 
4(>t> 

H07 

4rt; 

4fl9 

HU; 

HU 

H12 IP 

413 

HJ4 

415 i 7 

Hit, 

417 1 
Hit, 

Hl9 * 

H<ir> 

421 5 

422 
H,.>.J 4 
434 
H23 
H2fc> 

427 r 
426 r 
4 5; 4 r 
430 c 
H3i r 

432 C 

433 
4 34 
43s 
439 
437 
43b 

439 29 

440 


IF .0F-*W)t»O TO 1 

FNf;it'“V=2*r,|t * «H(NK » r. K | )*♦?/( ?»I 1< 4,2*IH1 ) 
f l\|OUi-’ = F filjt;. + fi OPF P 
UK*>' A ti = 

IK (ft.(! .1 , fit+31 ) 00 T 0 * 

fk =(I*' + ._**i + A. )^i14*X.)*,T\ + 2 *M+2. )t f-6 

COMTIMUK. 

Pl=C-l )**,(. f -l i*2»»;"*FAC*C**f*u(hK ,1 ) / ( (?*M+1 ) 

)**»t '•-! )M2*'»-l>*F„C*r**(K-l>/(2*FAC2>*3.l4159 
l*tFFAf.**2/sU)*«>*r /LUi r, f rjM 
Ph vsAUS(P >| 

IKfKR 2 .Pt , i. OF* 30) r,0 TO 4 
FK r.C= (■ '+■• +1 > »f pAC/ r \|-ri+? 1 
IF (FF/.r.b) .l.oF + l/,NAHAX = MK 
F 4 s r( 1 “ I X * K 2 
Aa = Al'3 ( FntihP J 
t“l = ALOPJ n ( M ) 

AH-A‘iS(l<2) 

A2 = A1.t.‘>l>i( AH) 

A3-A1+32 

IF «a3.9T.3n. )(.0 TO 3 
A3=/it'r.(i>i) 

A 1r/L ON)fl(A 5) 

AarA96<Ai) 4 Afi<;<A2) 

IK . r . .Mi . )t,0 TO 3 

PH ( f. J ;i) = t/(FM0R«*K4) 

( OfTIfj'4- 

K = t'+2* f -1 ) **HM*FH 1 t>,r;)*Ss0(rdX)*50K | 

COr (ij?j)=2»i> *K1*F4 (NM)+S|Jp. 

cO'irii'ji.'t 

PfrulMt 

FNnH.'lF'sO 
00 to a 

KH,ijmi = (0.,0.) 
t-u TO 16 
GK = n 

r 'U TO Id 
KOMaXbW 
14 = 1 
KtTUIO'l 
ENn 

f’U iVATIVtS OF the t<7(‘lAl FLt M C T TONS pi AND P4 
OF Oflt'ES K=n,M, A t.OIM Or APPL'meHT 0, 
also iv'ok^ali/atiom FurifTro^* n. 

Acpo C.JEKFK IFl'iT CnF FOP cOFPUIatIOM OF FO, 

SU'.ROOTlfJf. C [ID A OP < C «f’l t A • 10 t “HHaX • NNFIAX < tRRK A y) 

COrMlIv F IO( Hit)) »*’11 Of)»‘lO ) ,L ^rG)60) »F4 (59) t F4p(40) *CoF | 50) 
1 iSr ( ko » ,sof c-.o ) 

CU.PLLX 1X,PU,f4.RUH, F 4P, a , co f 
IX-(0.0,1,0) 

FFac=1 

A=(0.0.0.C) 

Ot) 17 NN=1,NNMAX 






44, l-ih = ?*rjN 

‘♦‘4,* ('■ = ?*Nf,,-l 

'♦ <+ SU"=0 

4'»4 f Nnnnsil, 

'445 00 If .'!|< = 1 . 1KRMAX 

T K =«? * i-J -1 

447 SO*- — '-j * t-V i, (t K « WK ) 

44.. ) 

449 f.\.OC--’U)G(, TO 1 

4? tr IY6.M4 -J-» ( L> ( ilk »W- ) ) *•<!/ t 2*JR+1 ) 

4b). ? t 14 Ml' 

45P ,6 COil| INI.'r 

4ji Hr = (-l >*• (' r-,_1 >*C«n(t!K,l )/(3*Sur ) 

4b‘4 l'ti>r < -1 » ** (• 0-1 I *3*4, 1'4 J «>«,/<r*C*0(NK.l) * S U M * 2 > * F F A. C 

45b M<^sAl»S(Ki^) 

4bf. IKjMIv.Ni ,i.(jr+30)f.O TO 4 

4 37 FK4C=C‘ + .')*44-ar/(0 + l) 

4‘jj; J F I 4-F (>,<-. ,4’t.l .of ■*•17 ) NFjMAXsMfc 

4b‘i f"»4' = Ki , *-j Xt'VP 

44.1. A&sftUS <FHCRM 

4t1 Al:AL:j r -lO( Ab) 

4h,' AHrAt'SIPf.P) 

I6v U ( A u ) 

464 A ir:A 1 4A2 

44-3 14 / /M. r -r,.A (J . )GD TO S 

404 Ab-AHS<Fi.4') 

44,7 Ai = AL'J u l»M A^) 

4ol- A 1= A^S ( A\) ♦«(T. ( AiJ ) 

U V 1 If ( Ai.f-T.5t, . >4,0 TO 3 

474 t4i>t r 

47, lK rU''tIIM'T 

4 7,' A )»*«r t ri-X >*F4PtflM»»SnP(NO)*Sl)l' 

47s cut (^,,>=cr(-l|*IX*ol t '*f4Pt( ; H)*S>Uf‘ 

47u ,/ (U.itiu'JL 

47b HtTlJ44N 

476 J PlvrirtHAin 

477 PU TO 2 

47J •« f 4(.,W|j ) = t 0 . ,0. >' 

47’ ftU TO If 

Hf,,i 4 

461 IU-1 

40;> PtTOKhJ 

44), fill) 

4ft, c 

4 j:, (- aScqC l.'.TtC 1-fGEA'C’PC FOI. Tf,tiM t US OF APOU^EIJT FTaO| 

4f*«, r uHOFT r-.rij AjTO i'" .1 iv/F’fi! P* UP to OPUf 4 w=P+P*NtJfi fll T"2l 

4fi/ c WIT. 4 W-!4 Ol/){ PO, tiP TO riK'.FP i'4aM+ii*T .'ir < lAX.l, 

4 (i r 

4,-9 SUl.kOUT iKfc Put YN(tTAO,f‘.Nf,MAX) 

40,, 01 E'-JSION P‘ J <3> 

491 CO.-woi-j 4-lf-t 1 00) tO, 10U.50, ,LfFGt50» «F4 (5o» ,F«i>(SO) .COF (50) 

4 1 .bi , ;,o ) ,j u, (f,fj) »Pt 50 ) ,PO(bO ) 

49 a rUf'PLh < F 4 »F 4p*CoF 

454 ?« = S )'< r ( 1-4 T AO^tTAr, , 

495 PP(T)sO 
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**96 PKtS>>3l 

4‘)/ Ih jO X 

49i >'0 2 U = li^ 

H-J j ? Pp I ) *sQ*Pp (2 ) 

500 1 Ml ) = PP ( 2 ) 

bni MU'’Axl = f Jl'-MAx + l 

bf :2 CO 3 lA’ = 2 ,Nr OaKi 

50 3 fsi=«- + J» ^ TJ N — a 

bo 1 . PO || t. -1« 2 

505 PP(b)=C (2*|-.-l )*E.T CO*PP{? > - ( N*«-l )*PP( 1) ) /( N-V) 

50S If (Wl.fD.'iNfnXH 00 TO 5 

507 

bOO PP(l)=l‘Ptii> 

b«9 4 PH(p> = i P<b) 

bl(> l'(iW>=PP(3) 

bir c COnTIi'jOE 
bl 2 Nl=IJN-1 

bib 3 PO(ML )-PHI 1) 

514 RETURN 

bib ewn 

blfc c 

bl 7 C. nui ATI. <;PHf PfJTPAL AivOuLAR K.’MCTrrf’S. £» OF AAfiufcENfs 

bio C C AMD MAiiS f-ftnEH r,\. With r-tf »-V£M; Sr t a • HP To ORrjER 

5is r w=-'+»*'iijNf.x-2« wi»h f'-> o'jo: s l taou* op to oropR 

b?u r tSfUiSH^UrtAX-l. 

52x c 

522 SOoKOUTlfJe Ofc aN(5 I WWMaX • IRH^iaX > 

523 CUK'MoiJ C H < 1 on ) «t.’C 1 UO »* 0 j ,unrG(5r > • ru ( bp ) , F4 P (sP ) , C;)F ( 50 ) 

524 l.bn(00),iCF(b0».K(no J,pn,bO),S lTA( bn)tSETAOt jbft) 

5^5 COvPLtX F4.F4p«C0F 

bpb N*n 

527 00 x i|N=l ,MjPaX 

520 NK-?*NM-X 

524 r-LTA('v M )S0 

5.Ju PK0s2*f-fJ 

531 5tTA(lu(;jO»=0 

532 PO 2 iM(=l,lKRnAX 

533 StTAOu(Mo) = ! 'tTAOr.i(M’j)+f-(AiA"«TPH)*PO(T K R) 

534 Scta< «\l> 1 J = SE Ti\( NN) ♦(;( wp,. Irh ) *PJ iKR » 

bob ’j COrrTl.jUK 

53b i COmTIUUC 

537 Rt-TURN 

533 ENn 

539 f 

540 C SUbROOTlNL HELD 0 

541 c 

542 SOrROUTIAt f IFL'JO C A TIP, A X • FU > 

543 COa.mOiv fJ<- Clop J »0 • lot » 5(1 ) »IJMFG(5c) t FA C 5n) , FUp (sc ) .COFC50) 

544 1 • br ( SO ) «s0( (50) »PCSP) tPOCbO) «St.TA(bO) »SeTA0tl5nl 

54b tO'PLLX (-OMF.I).F4 iF4P 

54c t0s(U.0,U.0) 

547 no 1 I|A; = 1,A‘C|H*AX 

54f, e 0=f 04COF (A f> ) *^ETAoQ (NM) 

549 1 CONTINUE 

550 RtruR^ 




SUb^ouiint HHOl 



